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We present in detail a simple, exact solution of the anisotropic 2-channel Kondo (2CK) model 
at its Toulouse point. We reduce the model to a quadratic resonant-level model by generalizing the 
bosonization-refermionization approach of Emery and Kivelson to finite system size, but improve 
their method in two ways: firstly, we construct all boson fields and Klein factors explicitly in terms 
of the model's original fermion operators Ckcrj, and secondly we clarify explicitly how the Klein 
factors needed when refermionizing act on the original Fock space. This enables us to explicitly 
follow the adiabatic evolution of the 2CK model's free-fermion states to its exact eigenstates, found 
by simply diagonalizing the resonant-level model for arbitrary magnetic fields and spin-flip coupling 
strengths. In this way we obtain an analytic description of the cross-over from the free to the 
non-Fermi-liquid fixed point. At the latter, it is remarkably simple to recover the conformal field 
theory results for the finite-size spectrum (implying a direct proof of Affleck and Ludwig's fusion 
hypothesis). By analyzing the finite-size spectrum, we directly obtain the operator content of the 
2CK fixed point and the dimension of various relevant and irrelevant perturbations. Our method 
goes beyond previous conformal field theory results, since it works for arbitrary magnetic fields 
and can easily be generalized to include various symmetry-breaking perturbations, and to study the 
crossover to other fixed points produced by these. Furthermore it establishes instructive connections 
between different renormalization group schemes such as poor man's scaling, Anderson- Yuval type 
scaling, the numerical renormalization group and finite-size scaling. 



I. INTRODUCTION 

Quantum impurity models displaying non-Fermi- 
liquid (NFL) behavior have attracted substantial in- 
terest during the past few years. These models 
have the common property that their exact elemen- 
tary excitations are not free-electron like and cannot 
be described using Fermi-liquid theory. Such single- 
impurity models have been proposed aSrjelevant for cer- 
tain properties oi,-heavy-fermion alloysocl and high-T c 
superconductors.Qu They aJse. emerge in the tunneling- 
impurity Kondo problem,EnlJ and infinite-dimensional 
strongly-correlated Lattice models can be mapped onto 
such models as well.liil All these models possess regimes 
in which the physical quantities are non-analytic (log- 
arithmic or power-law) functions of parameters such as 
temperature or magnetic field. 

The two-channel Kondo (2CK) model, introduced in 
1980 by Nozieres and Blandin,E3 is one of the simplest 
and most-studied quantum impurity models with NFL 
behavior. In this model two channels of spinful con- 
duction electrons interact with a single spin 1/2 impu- 
rity via a local antiferromagnetic exchange interaction. 
In contrast to the single-channel Kondo (1CK) model, 
which has a stable infinite-coupling fixed point at which 
the conduction electrons screen the impurity spin com- 
pletely by forming a spin complex, in the two-channel 
case they overscreen the impurity spin at infinite cou- 
pling, leaving a non-trivial residual spin object, so that 



the 2CK model's infinite-coupling fixed point becomes 
unstable. A stable fixed point exists at intermediate cou- 
pling strength, which leads to the appearance of a non- 
zero residual entropy and to non-analytical behavior for 
various physical quantities. Such non-analytical behavior 
was directly observed, for example, in anomalous conduc- 
tance signals of metallic-aiaaoconstrictions containing 2- 
state tunneling systems,L31lJ which are perhaps the most 
convincing realizations of 2CK physics found experimen- 
tally. 

The two-channel Kondo model has been stud- 
ied theoretically by an impressive number of differ- 
ent methods, which are comprehensively reviewed in 
Rcf. [l^. These m^ude-, approximate methods 
as the multiplicative3'E3o and the path-integral 
renormalization group approaches and slave-boson 
methods;EHrE-3 effective models such as the so-called com- 
pactified modeled eZI which is partially equivalent to 
the 2CK-, model; the numerical renormalization group 
(NRG);^~^2 and exact methods, such as_ .the JBcthe 
Ansatz,ErEl conforniaLfield theory (CFT)BEm and 
abelian bosonization.ae3~Ej Excepting abelian bosoniza- 
tion, however, the price for using powerful numerical or 
exact methods has hitherto always been a very high de- 
gree of technical sophistication and a lack of physical 
transparency. The Bethe Ansatz provides an analytical 
solution of the model, allowing for the calculation of the 
cross-over from Fermi-liquid to non-Fermi-liquid behav- 
ior of the thcrmodynamical quantities. It is, however, 
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rather involved, and is unable to calculate dynamical cor- 
relation functions. With the numerical renormalization 
group technique, which likewise is able to describe cross- 
over behavior, one can obtain thermodynamical proper- 
ties, carry out a finite-size analysis of the model, inves- 
tigate the effect of various perturbations such as differ- 
ent electronic and impurity magnetic fields and channel 
anisotropy, and in principle also calculate dynamical lo- 
cal correlation functions of the impurity. However, this 
method, though powerful, is approximate by construc- 
tion, requires considerable numerical prowess, is physi- 
cally not very transparent and is not well-suited to cal- 
culate dynamical properties of the conduction electrons. 
Finally, the elegant conformal field theory solution of Af- 
fleck and Ludwig (AL) focuses exclusively on the NFL 
regime in the vicinity of T = fixed point. By ex- 
ploiting its symmetry properties to the full, it provides 
the finite-size spectrum of the model, all thermodynami- 
cal quantities and furthermore all dynamical correlation 
functions. However, the CFT solution relies crucially on 
the so-called fusion hypothesis that can only be verified 
a posterio by comparing the CFT results with other ex- 
act methods. Moreover, it cannot be used to calculate 
the cross-over behavior, and requires, of course, extensive 
knowledge of the technical subject of boundary conformal 
field theory. 

A major advance towards finding a simple and trans- 
parent exact solution of the 2CK model was achieved 
by Emery and Kivelson (EK)ta with the rather simple 
technique of 1-dimensional abelian bosonization (peda- 
gogically reviewed in Ref. Q). Using bosonization and 
refermionization, EK showed that along a certain line 
in parameter space, known as the Toulouse "point" or 
Emery-Kivelson line, the anisotropical 2CK model can 
be mapped exactly onto a quadratic resonant-level model, 
which can be solved straightforwardly by diagonalization. 
Since spin anisotropy is known to be irrelevant for the 
multichannel Kondo model, this yielded new insight also 
about the generic behavior of the isotropic 2CK model. 
Though their approach works only in the vicinity of the 
EK line, the latter connects the Fermi-liquid and non- 
Fermi-liquid regimes, so that EK's method captures both 
the model's NFL behavior and the cross-over from the 
free to the NFL fixed point. EK calculated a number of 
thermodynamic and impurity properties and some elec- 
tron correlation functions, and related the NFL behavior 
to the fact that, remarkably, only "one half" of the im- 
purity's degrees of freedom (a Majorana fermion instead 
of a proper complex fermion) couple to the electrons. 

In the present workynwhich is an extended version of a 
previous publication,^ we generalize EK's bosonization 
technique to finite system sizes. For this purpose two 
important modifications are needed: 
(i) While EK use the field-theoretical approach to 
bosonization in which the bosonization relation ip a j ~ 



F, 



-i(fca 



the boson fields <fi a j and Klein factors F a j are constructed 
explicitly from the original ip a j operators, so that the 
bosonization formula becomes an operator identity in 
Fock space. 

(ii) Since EK were interested mainly in impurity proper- 
ties, they did not need to discuss at all the Klein factors 
F a j [which lower the number of aj-electrons by one and 
ensure proper anticommutation relations for the ij) a j's]. 
These Klein factors, however, are essential for quantities 
like the finite-size spectrum or various electron correla- 
tion functions. Therefore it is crucial to specify how the 
Klein factors for the refermionized operators act on the 
Fock space. As we shall see, these new Klein factors are 
only well defined on a suitably enlarged Fock space that 
also contains unphysical states, which must be discarded 
at the end using certain gluing conditions. 

With these modifications, EK's bosonization approach 
enables us by straightforward diagonalization of the 
quadratic resonant-level model (i) for the first time to an- 
alytically trace the cross-over of the 2CK model's finite- 
size spectrum from the FL to the NFL fixed point, at 
which we reproduce the fixed-point spectrum previously 
found by CFT using a certain fusion hypothesis; (ii) to 
construct the eigenstates of the 2CK model correspond- 
ing to this crossover spectrum explicitly; (iii) and to ex- 
tract the operator content of the NFL fixed point and 
determine the dimensions of different relevant and irrele- 
vant operators. Since our method works also in the pres- 
ence of an arbitrary magnetic field (unlike CFT), we can 
also (iv) investigate how a finite magnetic field destroys 
the NFL spectrum for the low-energy excitations of the 
model and restores the FL properties, (v) Furthermore, 
our finite-size bosonization approach can easily be re- 
lated to various popular renormalization group methods; 
it therefore not only provides a useful bridge between 
them, but can potentially be used as a pedagogical tool 
for analytically illustrating-their main ideas. 

In a future publication^ we shall show that this 
method furthermore allows one (vi) to construct very 
easily the scattering states of the model; (vii) to prove 
explicitly the validity of the bosonic description of the 
NFL fixed point Maldacena and Ludwigj!3 (viii) to de- 
termine the fixed point boundary conditions at the im- 
purity site for the different currents and fields in a very 
straightforward way; and (ix) to calculate with ease all 
correlation functions at and around the NFL fixed point. 
This implies that all CFT results can be derived from 
first principles using the bosonization approach. 

The paper is organized as follows. In Section [l| we de- 
fine the 2CK model to be studied. For completeness, 
and since the proper use of Klein factors is essential, 
Section III briefly reviews the "constructive" (operator 
identity-based) approach to finite-size bosonization used 
throughout this paper. The Emery-Kivelson mapping 
onto a resonant-level model is discussed in Section IV, 



is used merely as a formal correspondence, we 



use the more careful and , 
tion procedure of Haldand 



cplicit constructive bosoniza- 
I. In the latter approach both 



using our novel, more explicit formulation of refermion- 
ization within a suitably extended Fock space. The solu- 
tion of the resonant level model and the construction of 
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generalized gluing 
In Section VI the 



the non-Fermi liquid spectrum usin 
conditions is presented in Section 
results of our finite-size calculations are compared with 
and interpreted in te rms of various RG procedures. Fi- 
nally, in Section VII we summarize our conclusions and 
compare our method to a few others. 

The centerpiece of our work, and indeed the prereq- 
uisite for all of our results, is the uncommonly care- 
ful and detailed finite-size formulation of the Emery- 
Kivelson mapping in the main text. Technicalities not 
related to this mapping are relegated to the Appendices: 
Appendix jX] discusses in some detail matters related to 
the choice of an ultraviolet cutoff, and also gives the 
often-used position-space definition of the 2CK model, 
to facilitate comparison with our momentum-space ver- 
sion. The construction of the extended Fock space needed 
for refermionization is discussed in Appendix [b|, and the 
technical details used to diagonalize the resonant-level 
model and to calculate several of its properties are given 
in Appendix ^|. Finally, for pedagogical reasons and for 
the sake of completeness, in Appendix [d] we use our 
finite-size bosonization method to solve the 1-channel 
Kondo model as well. 



II. DEFINITION OF THE MODEL 

A. Hamiltonian in Momentum Space 

Throughout the main part of this paper we shall use 
the standard 2CK Hamiltonian in momentum space (its 
position-space representation is given in Appendix 0). 
We consider a magnetic impurity with spin 1/2 placed at 
the origin of a sphere of radius R — L/2, filled with two 
species of free, spinful conduction electrons, labeled by a 
spin index a = (T,J.) = (+, — ) and a channel or flavor 
index j = (1,2) = (+, — ). We assume that the interac- 
tion between the impurity and the conduction electron is 
sufficiently short-ranged that it involves only s-wave con- 
duction electrons, whose kinetic energy can be written as 
(v F = H = 1) 



The operator cj_ aj - creates an s-wave conduction electron 
of species (aj) with radial momentum k = p — pf rela- 
tive to the Fermi momentum pp, and the dispersion has 
been linearized around the Fermi energy sf- Sk — £F ~ k. 
The symbol : : in Eq. ([j]) denotes normal ordering with 
respect to the free Fermi sea or "vacuum state" |0) (the 
reason for this notation will become clear in the next 
section), defined by 



J2 k 

kaj 



: C 



kaj C kaj ■ i 



(1) 



Cfcaj|0) =0 for k > 0, 
ct aj \0) o = for fc<0. 



(2a) 
(2b) 



The Cfccy 's obey standard anticommutation relations, 

{c&etj ^U'cc'i' } — &kk' $aa'&jj' 



33 ' 



(3) 



where due to radial momentum quantization in the spher- 
ical box, the values taken on by k are quantized: 



-(n fe -Po/2), n k eZ 



(4) 



Here Po = or 1, since at zero temperature the chemi- 
cal potential (and hence pf) must either coincide with a 
degenerate level (Po = 0) or lie midway between two of 
them, respectively (Po = !)■ Evidently the level spacing 
in both cases is 



= T' 



(5) 



Since the s-wave conduction electrons form an effec- 
tively one-dimensional system they can alsoJae described 
by a one dimensional chiral field, defined asE3 



E 



{ip a j(x),ip a >f(x')} = 6 aa >S jjr 27T5(x - x') . (7) 

In the continuum limit L — > oo, the x > and x < por- 
tions of (x) can be associated with the incoming and 
outgoing scattering states, respectively. Note that for 
Po = or 1 the fields ?p a j(x) have periodic or antipcriodic 
boundary conditions at x — ±L/2, respectively, hence Po 
will be called the "periodicity parameter" . 

We assume a short-ranged anisotropic exchange inter- 
action between the impurity spin and the s-wave con- 
duction electron spin density at the origin, i.e. a Kondo 
interaction of the form 



Pir 



[i,k,k' 



A^^t : c fcQ , J -( 2 cr„ Q /)cfe' Q ,' J ■: 



(8) 



Here the — x, y, z) are the impurity spin opera- 

tors, with S z eigenvalues (t|",JJ-) = an d the A M 's 

denote dimensionless couplings: A 2 generates different 
phase shifts for spin-up and spin-down conduction elec- 
trons, while \ x = X y = X± describe spin-flip scattering 
off the impurity. Finally, the effect of a finite magnetic 
field is described by 



H h = hS z + h e Ns, 



(9) 



where hi and h e denote the magnetic fields acting on the 
impurity and conduction electron spins, respectively, and 
J\ s (to be defined slightly below) denotes the total spin 
of the conduction electrons. 

Finally, note that we have taken all sums over 
fcrmion momenta above to be unbounded, since the con- 
structive bosonization scheme we intend to use requires 
an unbounded fermion momentum spectrum. We thus 
have effectively taken the fermion bandwidth, say D, to 
be infinite, but will reintroduce an ultraviolet cutoff when 



defining the boson fields in Eq. (14) below 
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III. BOSONIZATION BASICS 

The key to diagonalizing the Hamiltonian is to find 
the relevant quantum numbers of the problem and to 
bosonize the Hamiltonian carefully. While the tech- 
nique of bosonization is widely used in the literature, 
the so-called Klein factors mentioned in the introduc- 
tion are often neglected or not treated with sufficient 
care. However, it, has recently been emphasized by sev- 
eral authorscJ'OL30 that these Klein factors are very 
important in some situations, one of which is the calcu- 
lation of the finite-size spectrum. In the present Section 
we therefore discuss our bosonization approach in some- 
what more detail than usual, formulating it as a set of 
operator identities in Fock space, and emphasizing in par- 
ticular the proper use of Klein factors to ladder between 
states with different particle numbers in Fock space. (An 
elementary, pedagogical and detailed introduction to the 
bosonization scheme used here, which is based on that of 
Haldane,c3 may be found in Ref. ||.) 



A. Bosonization Ingredients 

To characterize the electronic states, we start by intro- 
ducing the number operators 



C kaj Cka i ' ' 



(10) 



which count the number of electrons in channel (aj) with 
respect to the free electron reference ground state |0)o- 
The non-unique eigenstates of N a j will generically be de- 
noted by \N) ee \N n ) ® \N±i) (gi |A T2 ) ® |iV i2 ), where the 
N a j's can be arbitrary integers, i.e. N G Z 4 . 

Next, we define bosonic electron-hole creators by 



fct 



c k+ qa j c k<xj , (q = 2nn q /L > 0) 



n k eZ 



(11) 



<ia.) 



where the n q are positive integers. The operators b 
create "density excitations" with momentum q in chan- 
nel aj, satisfy standard bosonic commutation relations, 
and commute with the N a j's: 



l b qaj' b l> a >j>} - S qq' S c, a >Sjj> , 

[b qaj ,N a , f ]=Q. 



(12) 



Among all states \N) with given N, there is a unique 
state, to be denoted by |7V) , that contains no holes and 
thus has the defining property 



b qa j\N) = (for any q > 0, a,j) 



(13) 



We shall call it the 'W-particle ground state" , since in the 
absence of the interaction term (ph, no \N) has a lower 



energy than \N) ; likewise, no \N)o nas a lower energy 
than the "vacuum state" |0)q defined in Eq. (0). Note, 
though, that if P = 0, the states co aj |0)o are degenerate 
with 1 0) o , because then co a j removes a zero-energy elec- 
tron. [In the conformal field theory literature the states 
\N)o are sometimes referred to as [J7(l)] 4 primary states, 
since the quantum numbers N a j are just the charges 
associated with the ?7(1) gauge transformations of the 
fields tp a j — ► tfj aj e l&ai .] It can be provencalfil that any 
TV-electron state \N) can be written as \N) = fQy)\N)o, 
i.e. by acting on the A^-electron ground state with an 
appropriate function of electron-hole operators. 
Next, we define bosonic fields by 



cj) aj (x) = Y 

q>0 



-1 



JqxiA 
- u qaj 



(14) 



Here a ~ 1/pf is a short-distance cutoff; it is intro- 
duced to cure any ultraviolet divergences the theory may 
have acquired by taking the fermion bandwidth D to 
be infinite. It is well-known, however, that within this 
bosonization cutoff scheme the coupling constants have 
different meanings than for other standard regularization 
schemes using a finite fermion bandwidth, and that the 
relations between coupling constants in different regular- 
ization schemes can be found by requiring that they yield 
the same phase shifts. For the sake of completeness, we 
discuss this and other cutoff related matters in some de- 
tail in Appendix 

It is easy to prove that the fields d x 4> a j{x) are canoni- 
cally conjugate to the <p a j(xYs, in that 



[<f) aj (x),d x '(p a 'j'(x')] = 2iri(6 a (x- x') - 1/L) 5 aa >Sj 

where S a (x — x') is a smeared delta function: 

a/ir 



(15) 



S a (x - x') 



{x — x') 2 + a 2 



(16) 



As final bosonization ingredient, we need the so-called 
Klein factors F a j, which ladder between states with dif- 
ferent Naj's (which no function containing only I r qa - can 
accomplish, since these conserve N a j). By definition, the 
Faj's are required to satisfy the following relations: 

[■F 'xj 3 N a fjf ] &a.a'^jj'F a j , (17a) 

[F aj ,b qa/j ,] = lF aj ,b qalj ,] = 0, (17b) 

V<ti = F ij F *j ^ 1 - (17c) 
{F a] ,Fl J ,} = 2 5 aa ,5 jr (17d) 
{F aj ,F a , f } = for (aj) ± (a'f) . (17e) 

These relations imply that when F a j is applied to a state 



\N) = /(&t)|jV) , it commutes past /(6+) [by (17b)], and 
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then removes [by fll7a| )] an (aj) electron from the top- 
most filled level of |7V)o, namely Uk,aj — N a j\ to be ex- 
plicit, F a j\N) = f(b^)cN a . a j\N) . Thus F a j decreases 
the electron number in channel {aj} by one, N a j — > 
N a j — 1, without creating particle-hole excitations. Sim- 
ilarly, F^j adds a single (aj) electron. As shown in 

Refs. ||or ||, the construction F aj = a 1 / 2 tp aj (0)e l4 '^ < - \ 
which explicitly expresses F a j in terms of the fermion op- 
erators Ckaj, has all the desired properties. 



IV. MAPPING ONTO RESONANT-LEVEL 
MODEL 

In this section we map the 2CK model onto a reso- 
nant level model, using a finite-size version of the strat- 
egy invented by Emery and Kivelson: using bosonization 
and refermionization, we make a unitary transformation 
to a more convenient basis, in which the Hamiltonian is 
quadratic for a certain choice of parameters. 



B. Bosonization Identities 



Having introduced the Klein factors F a j and the bo- 



son fields 



we are ready to bosonize, i.e. to rewrite 



expressions involving the fermion operators Ck a j in terms 
of the bosonization ingredients defined abovfe-. In our 
notation, the standard bosonization identities^ for the 
fermion field, density and kinetic energy take the follow- 
ing forms :C3 

lf>aj(x) = Fa]a -^ e ~^ a3 -P /2)2, x /L e -i^(x) ^ (18) 

:^ aj {x)^ aj (x) := ^-d x< f>aj(x) + N aj /L , (19) 

Z7T J Z7T 

Ho = ^ aj (N aj + 1-P ) + J2q b\ aj b qaj . (20) 

aj aj 

q>0 



Several comments are in order: Firstly, in the limit 
a — > Eqs. ( |l8| ) to ( po|) are not mere formal correspon- 
dences between the fermionic and bosonic expressions, 
but hold as rigorous operator identities in Fock space. 
For a ^ 0, ( p^|) and (|l^) are not rigorously exact, but in- 
stead should be viewed as conveniently regularized redef- 
initions of th e fer mion fields and densities (as discussed 
in Appendix 



factors F a j p. 



A 2). Next, in Eq. (g 



for ip a j , 



the Klein 



ay a twofold role: firstly, by Eq. (17a) they 
ensure that the right-hand side of Eq. ( |l8| ) acting on any 
state indeed does lower the number of aj-electrons by 
one, just as ip a j does; and secondly, by Eqs. (17d) and 
(17c) they ensure that fields with different (aj)'s do have 
the proper anticommutation relations (Q). In contrast, 
Eq. (|l^) for the den sity operator contains no Klein fac- 
tors [because of (17c)]. Finally, in Eqs. ( po[ ) for the kinetic 
energy, the first Al term is just o(N\H n \ N)o, the energy 
of the A^-particle ground state \N)o relative to |0)o- Since 
the Klein factors do not commute with this term, they 
evidently cannot be neglected when calculating the full 
model's finite-size spectrum, for which all terms of order 
Al must be retained. The second term of ( p0| ) describes 
the energy of electron- hole excitations relative to \0)o- 
[Its form can be obtained by observing that the commu- 
tator [bq a j, Ho] is the same when calculated in terms of 
Cfcaj's using 
and (|2T 



lh and (ILlh, or in terms of bq^s using 1l 



A. Conserved Quantum Numbers 



The quantum numbers N a j of Eq. ([L0|) are conserved 
under the action of Hq, and H z (the A 2 term of 
Hi n t = H z + H±), but fluctuate under the action of the 
spin- flip interaction H± (the A_l term). On the other 
hand, the total charge and flavor of the conduction elec- 
trons is obviously conserved by all terms in the Hamilto- 
nian, including H±. Therefore it is natural to introduce 
the following new quantum numbers: 



(M c \ 




f 1 


1 


1 


M s 






-1 


1 


M f 


■il 




1 


-1 


\mJ 






-1 


-1 




(21) 



where 2JV C , M s , and Mf denote the total charge, spin, 
and flavor of the conduction electrons, and M x measures 
the spin difference between channels 1 and 2. Clearly, 
any conduction electron state \N) can equally well be 
labeled by the corresponding quantum numbers M = 
(Mc,M s ,Mf,M x )- However, whereas the A^-'s take arbi- 
trary independent integer values, the M's generated by 
Eq. ( pl| ) (with iV e Z 4 ) can easily be shown to satisfy 
the following two constraints, to be called the free gluing 
conditions: 



M e [1 + P/2) 4 , 

M c ± Mf = (Ms ± A4) mod 2 



(22a) 
(22b) 



where the parity index P equals or 1 if the tot al nu mber 
of electrons is even or odd, respectively. Eq. (22a) for- 
malizes the obvious fact that the addition or removal of 
one aj electron to or from the system necessarily changes 
each of the M y 's by ±1/2, so t hat t hey are either all in- 
tegers or all half-integers. Eq. (221) selects out from the 



set of M of the form (22a) those that correspond in the 
old basis to TV e Z 4 , i.e. to physical states (an M of 
the form ( p2a| ) that violates ( (22b| ) would correspond to 
G (Z + l/2) 4 , which does not exist in the physical Fock 
space) . 

The new basis has two major advantages: firstly, M c 
and Mf are conserved quantum numbers; and secondly, 
the quantum number M s fluctuates only "mildly" be- 
tween the values St T 1/2, since the total spin, 
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iK + s, 



(23) 



is conserved. In contrast, the quantum number M x fluctu- 
ates "wildly", because an appropriate suc cessi on of spin- 
flips can produce any J\f x that satisfies ( |22b| ), as illus- 
trated in Fig. [l| This wildly fluctuating quantum number 
will be seen below to be at the heart of the 2CK model's 
NFL behavior. In revealing contrast, the 1CK model, 
which shows no NFL behavior, lacks such a wildly fluc- 
tuating quantum number (see Appendix . 

Since St, M c and TV/ are conserved, the Fock space 
•^phys of all physical states can evidently be divided as 
follows into subspaces invariant under the action of H: 



phys 



@'S T ,M C ,M S 



S phys (S T ,Mc,Mf) 



(24) 



S phys (S T ,Af c ,Mf) = {Wc,S T -i/2,M f M x -A) 

®\K,S T + l/2,JVf,M x + l-A)}. (25) 



In both equations the prime on the sum indicates a re- 
striction to those Afy's that satisfy the free gluing condi- 
tions (|22]). To diagonalize the Hamiltonian for given St, 
Af c and Aft, it evidently suffices to restrict one's attention 
to the corresponding subspace S p h ys (ST,Af c ,Aff). 



B. Emery-Kivelson transformation 

Following E mery and Kivelson, we now introduce, in 
analogy to Eq. (pi|), new electron-hole operators and bo- 
son fields via the transformations, 



<py — iZocj Ry,aj4>aj 



(y = c,s,f,x) , (26) 



where Ry, a j is the unitary matrix in (21). These obey 



relations analogous to (|T^) and (|l5|), with aj — > y. More- 
over, we define \Af)o, the A/"-particle vacuum state, to 
satisfy b qy \N)o — 0, as in (|l3|). If M and N are related 
by (|2l|), then the states \A7)o an d \N)q are equal up to 
an unimportant phase (see Appendix ^), because both 
have the same N a j and Af y eigenvalues and both are an- 
nihilated by all bqaj's and b qy 's. 

Using the quantum numbers Af y and the bosonic fields 
tp y (x), the Hamiltonian takes a transparent form. The 
free electron part Hq of (^0|) can be written as 

H a = A L [aT c (1 - P ) + J2^y/ 2 } + E 1 h \v h iv . 

y y, <?>o 

(27) 

while Eqs. © and (||) are used to obtain, respectively, 



H z = X z [d x tp s (0) + AjX] S z , (28) 



H 



2a 



e -^ s (o) 5+ ( F t F ne -i<p 



(0) 



+ J F, t 2 F T2 e^ (0) )+h, 



(29) 



Eqs. ( |27j ) to ( |2S| ) and (^|) constitute the bosonized form 
of the Hamiltonian for the anisotropic 2CK model, up to 
and including terms of order At ■ 

Next we simplify H Zl which, being diagonal in spin 
indices, merely causes a phase shift in the spin sector. 
This phase shift can be obtained explicitly using a uni- 
tary transformation (due to EK) parameterized by a real 
number 7, whose value will be determined below: 



H -> H' = UHU^ 



(30) 
(31) 



Under this transformation the impurity spin, the spin- 
diagonal part of the Hamiltonian, the spin boson field 
and the fermion fields transform as follows (rising e.g. 
the operator identities of Appendix C of Ref . Q) : 

S± -> US±U^ = e ±<w »(°)5± , (32) 
H + H z -> H + (\ z - j)d x <p s (0)S z 

+ \ Z A L N S S z + 7 2 [l/(4a) - tt/(4£)] , (33) 

ifs(x) — > (fi s {x) — 2^S Z arctan(x/a) , (|.x| <C L) . (34) 

i; aj (x) -> Vai(^)e ia7Sa "cta„(*/a) , (|x| « L) . (35) 

Eq. ([53]) is most easily derived in the momentum-space 
representation [using ([l2|), ( |l4| ) and (p0|), see Section 7 
of Ref. 44 1; on the other hand, since in Eq. ( ]34"|) we only 
give the |x| <C £ limit in which order l/L terms can be 
neglected, this equation is easier to derive in the position- 
space representation (by first evaluating Ud x j£ s (x)U~ 1 
using (|l5| ) and (|l6|), then integrating). Eq. ( |35|) follows 
from (||), since ^joc e ~ laip ^ 2 [by @ and ([2§)]. 

Recalling [from Pjj)] that d x ip s (x)j2Tr contributes to 
the conduction electron spin density, we note by dif- 
ferentiating (HJ) that the EK transformation produces 
a change in the spin density of — 2jS z ■n8 a {x) /2tt] intu- 
itively speaking, it ties a spin of —"fS z from the conduc- 
tion band to the impurity spin S z at the origin. 

To eliminate the S z d x (p s term in ([53]), we now choose 
7 = X z ; then the spin- flip- independent part of the Hamil- 
ton takes the form 

H'{X X = 0) = \ z A L Af s S z + ^Af y /2 

v 

+ E q b ly b iy + Hh + const ' ( 36 ) 

y,i>o 

and H' ± contains the factors e ±*(l-A»)*>.(0)_ These fac- 
tors are simply equal to 1 at the special line X z — 1, the 
so-called Emery-Kivelson line, at which H' ± simplifies to 



G 



2a 



S+iFl^xe-^W + ft F T2 e < *'-( )) + h.c. 



(37) 



C. Refermionization 



Definition of New Klein Factors 



We shall henceforth focus on the case A z = 1, which will 
enable us to diagonalize the model exactly by refermion- 
ization. Deviat ions from the EK line will be shown 
in Section VI O] to be irrelevant, by taking 7=1 but 
A 2 = 1 + 5X Z , and doing perturbation theory in 



6H' Z = 5\ t [a x <p.(p) + A L M S }S Z 



(38) 



Before proceeding, it is instructive to interpret the sig- 
nificance of the EK line using the standard heuristic lan- 
guage of Nozieres and Blandin. They argued that in the 
strong- coupling limit, an anti-ferromagnetically-coupled 
impurity will "capture" one spin 1/2 from both the j = 1 
and 2 channels (capturing a spin from just one channel 
would break channel symmetry), i.e. it will tie a spin of 
— 2S Z from the conduction band to its own S z , which 
yields an "overscreened effective spin of — S z n . However, 
since the latter is again anti-ferromagnetically coupled to 
the remaining conduction electrons, this strong-coupling 
fixed point is unstable in the RG sense, just as the weak- 
coupling fixed point. Thus, they argued that there must 
be a NFL fixed point at intermediate coupling which is 
stable (for h e = hi = 0). 

Now, a crucial property of the EK line is that it con- 
tains this NFL intermediate-coupling fixed point. A 
heuristic way to see this it to note that on the EK line, 
the impurity spin is in fact "perfectly screened" : the spin 
— ^Sz from the conduction band [mentioned after (pi])], 
that is tied to the impurity by the EK transformation, is 
equal to — S z if 7 = X z = 1, thus precisely "canceling" 
the impurity's spin S z . Thus, on the EK line the impu- 
rity "captures" exactly one spin 1/2 from the conduction 
band to form a "perfectly screened singlet" with zero to- 
tal spin (which is the heuristic reason why the EK line is 
stable), but it does so without breaking channel symme- 
try, since d x ip s is constructed in equal amounts from all 
four species of 4> a j fields. 

Of course, there are also more rigorous ways of seeing 
that the NFL fixed point lies on the EK line: Firstly, 
for X z = 1 it follows from (35) that the phase shift S of 
the outgoing relative to the incoming fields, defined by 
VW(CT) = e i2S ip aj (0+) (with |0±| > a), is \S\ = tt/4, 
which is just the valuc-kspwn for the NFL fixed point 
from ot her ap proaches EJEIJ Secondly, we shall deduce in 
Section VI C from an analysis of the finite-size spectrum 
that the leading irrelevant operators (with dimensions 
1/2) vanish exclusively along this line, but not away from 
it. Since the presence or absence of the leading irrelevant 
operators strongly influences the low-temperature PfiWr 
erties of the model (such as its critical exponents)H'E2l, 
and since these must stay invariant under any RG trans- 
formation, one concludes that the Emery-Kivelson line 
must be be stable under RG transformations. 



The most nontrivial step in the solution of the model is 
the proper treatment of Klein factors when refermioniz- 
ing the transformed Hamiltonian. In their original treat- 
ment Emery and Kivelson did not discuss Klein factors 
at all and simply identified e~ t!px ^ x ' j\fa as a new pseud- 
ofermion field ip x (x). Although this procedure happened 
to be adequate for their purposes, the proper consider- 
ation of the Klein factors and gluing conditions is es- 
sential, as already emphasized in the introduction, for 
giving a rigorous solution of the problem and obtaining 
the finite-size spectrum. Some other authors tried to im- 
prove the Emery-Kivelson procedure by representing the 



Klein factors by F a 



-i&a 



, where <d a j is a "phase op- 



erator conjugate to N a j" , and added these to the bosonic 



fields cf> a j before making the linear transformation (|2^ 
This procedure is problematic, however, since then the 
e -i¥> y (o)> g con t am f ac tors such as e -40Qj / 2 , which are ill- 
defined (for a more detailed discussion of this point, see 
Appendix D.2 of Ref. g). 

A novel, rigorous way of dealing with Klein fac- 
tors when refermionizing was presented in Ref. (and 
adapted in Ref. |5J to treat an impurity in a Luttinger 
liquid): We introduce a new set of ladder operators Tl 
and T y (y = c, s, /, x) to raise or lower the new quantum 
numbers Af y by ±1, with, by definition, the following 
properties [in analogy to Eqs. jl7|)]: 



[j-y,J\fyl] — Syy'J-y , 

[Fy,b qy >] = [T y ,b\ yl ] = , 

t ■rt — art 77 _ 1 

y v v v ' 

{Ty, J~yl} — 2 Syy' , 

{^,^}=0 fory^y'. 



(39a) 
(39b) 
(39c) 
(39d) 
(39e) 



Now, note that the action of any one of the new Klein 
factors Ty or T^ respects the first of th e free gluing con- 
ditions, (22a), but not the second, (22b). More generally, 
( |22b| ) is respected only by products of an even number 
of new Klein factors, but violated by products of an odd 
number of them. This implies that the physical Fock 



space .T-phys of all \J\f) satisfying both (22a) and ([22b]) 
is closed under the action of even but not of odd prod- 
ucts of new Klein factors. For example, let \ip) p hy S be in 
■^"phys, then TlTy'\ijj) phys is too, but T y \ip} phys violates 
(E2q) and hence is an unphysical state. The action of 
arbitrary combinations of new Klein factors thus gener- 
ates an extended Fock space T e xt, which contains ^" p hys 
as a s ubsp ace and is spanned by the set of all \M) satis- 
fying (22a), including unphysical states violating (22b). 
To demonstrate that ^phys can indeed be so embedded 
in T cx t, we explicitly construct a set of basis states for 
the latter in Appendix |b| 
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Since odd products of TyS lead out of .Fphys, they 
cannot be expressed in terms of the original Klein fac- 
tors F a j, which leave .Fphys invariant. However, and 
this is crucial, the Hamiltonian contains only even prod- 
ucts of old Klein factors. Now, any "diagonal" combi- 
nation F^F a j — 1; and any "off-diagonal" combination 
F^jF a ij> or FajF^tji acting on any state \N) just changes 
two of its N a j quantum numbers. Using (|2l]) to read off 
the corresponding changes in J\f s , J\ff and J\T X , we can 
thus make the following identifications between pairs of 
the old and new Klein factors: 



• r x- r s — ^fl^il 



T X T\ = -^12-^12 



f 



fIf 



pt pt 



(40a) 
(40b) 



Each of these relations involves an arbitrary choice of 
phase, whose consequences for the basis states \M) are 
discussed in Appendix [b|. These choices uniquely fix 
the phases of all other similar bilinear relations be- 
tween old and new Klein factors, which can be found 
from composing (and conjugating) the above four, e.g. 
T\T) = -{T x Tt)(TlT)) = FjjF^. Since the relations 
([h]) by construction respect ( pi| ) (as can be verified by 
acting on any |iV)), they, and all similar bilinear relations 
derived from them, also respect both of the free gluing 
conditions (|22]). The relations involving T c and Tf are 
not needed for the present 2CK model, but are included 
for completeness; for example, Tf would be needed for 
models involving "flavor-flip" processes. 

We can thus replace the Kle in fa ctor pairs occurring 
in Eq. @ by the ones in Eq. (ff0a|): 



^ 2a 



S+F s (F x e- llp * (0) + Fjje^-^) + h.c. 



(41) 



The only consequence of this change is that we now work 
in the extended Fock space .F ex t, and will diagonalize H' 
not in the physical invariant subspace S V Y iys {ST,Mci-!^'f) 
of (|B|), but in the corresponding extended subspace 
^ext^Tj A/" c , A//), given by an equation similar to (25), 
but w here the (B 'Af x sum now is restricted only to satisfy 
fl22a[ ), not also (22b). At the end o f the calculation we 
shall then use the gluing condition ( |22b| ), satisfied only 
by the physical states in <S GXt but not by its additional un- 
physical states, to identify and discard the latter (see Sec- 
This approach is completely analogous to the 



VC) 



tion 

use of gluing conditions in AL's CFT solution of the 2CK 
model. It is also in sonm-sense analogous to Abrikosov's 
pseudofermion techniques; there an impurity-spin oper- 
ator is represented in terms of pseudofermion operators 
acting in an enlarged Hilbert space, which contains not 
only the physical one-pseudofermion states, but also un- 
physical many- or no-pseudofermion states that are pro- 
jected away at the end of the calculation. 



2. Pseudofermions and Refermionized Hamiltonian 

We now note that H'^ of (f4l| ) can be written in a form 
quadratic in fermionic variables, 

H ' ± = ^(^(0)+^(0))(c d -4), (42) 

by defining a local pseudofermion and a pseudofermion 
field ip x (x) by the following refermionization relations: 



F*S- , c\c A = S Z + 1/2 , 



(43) 



il> x {x) = i^ a -V2 e -i(^ -1/2)2^/1, e -W») ( 44a ) 

(44b) 



-ikx 
e c ka 



where Eq. ( |44b| ) defines the cj, x as Fourier coefficients of 
the field ip x (x). For reasons discussed below, the field 



ip x in Eq. (44a) has been defined in such a way that 



its boundary condition at ±L/2 is P-dependent, since 
JV X € Z + P/2 and <p x (x) is a periodic function. Thu s 
the quantized k momenta in the Fourier expansion ( |44b| ) 
must have the form 



k = A L [n- k -(l-P)/2] 



(n s G Z) 



(45) 



(i.e. the periodicity parameter Pq of (Q) here equals 

The new pseudofermions were constructed in such 
a way that they satisfy the following commuta- 
tion/anticommutation relations 



{ c k x > c \' x } — &kk' j 

{cd, 4} = 1 , 
{cd,c\ } = {c d ,c- kx } = , 



Cd , 



(46) 
(47) 
(48) 
(49) 



which follow directly from the prop erties of ip x [by anal- 
ogy to the relations of Section [n| and Eqs. (|39|). Note 
that Cd lowers the impurity spin, raises the total electron 
spin J\ s and hence conserves the total spin St, whereas 
ip x conserves each of the impurity, electron and total 
spins. 

Next we have to relate the number operator for the 
new rr-pseudofermions to the quantum number Af x . This 
requires defining a free reference ground state, say |0)s ext , 
in the extended subspace iS ex t ! with respect to which the 
number of pseudofermions are counted. In analogy to 
(§), we define |0) 5ext by 



Cfcj0)5 oxt =0 for fc>0, 

c Ll°>^ t =0 for fc<0, 

c d |0) 5cxt = for e d > 0, 

4|0)5 oxt =0 for £d <0, 



(0) 



4 0) - 



(50a) 
(50b) 

0, (50c) 
= 1. (50d) 
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Here ej, whose value will be derived below [see (pq)], is 

the energy associated with the Cd pseudofermion, and n d 
denotes its occupation number in the reference ground 
state |0)s cxt . Using the symbol : : to henceforth de- 
note normal ordering of the pseudofermions w.r.t. |0)s 6Xt , 

we have : ctc^ ■— c\cd — n d °^ . Furthermore, we de- 
fine the number operator for the x-pseudofcrmions by 

N x = J2k -c\ x c- kx :. Then equations @, © and © 
together imply that 



■ N x - P/2 



(51) 



holds as an operator identity. This can be seen intu- 
itively by noting that ip x ~ F x ~ c kx [by (0)], hence 
the application of ip x (or tp x ) to a state decreases (or 
increases) both M x and N x by one. These two numbers 
can thus differ only by a constant, which must ensure 
that N x is an integer. Our definition of |0)s oxt effectively 
fixes this constant to be P/2, by setting N. x = fps. 
Af x — P/2. (This can be verified rigorously by checkings 

that lim^o lim a ^ /f0 2 ^[i>U x + x o)i>x{x) ~ 
when evaluated using either (|44a| ) or ( |44b| ) , yields the 
right- or left-hand sides of d51|) , respectively. Sim- 
ilarly, Eq. ( |5^ ) below can be proven by evaluating 

lim xo ^ lim Q ^ /f// 2 ^[^U x + x Q )id x ip x (x) - ^].) 

We are now ready to refcrmionize the Hamiltonian W . 
The kinetic energy of the k pseudofermions obeys the 
following operator identity: 



k :c\„c ka 



A, 



k 



^kx 



N X (N X + P) + Y,qbl x b qx . (52) 



This follows by analogy with (Q) and (|2(]), with N x 
and 1 — P instead of N a j and Pq. Now note that 
N X (N X + P) = A/; 2 - P/4, i.e. (||) does not contain 
a term linear in N x . Actually, the choice of the phase 



-i{N x -l/2)2irx/L 



in our refermionization Ansatz (44a) for 



ip x (x) was made specifically to achieve this. Hence ( |52j ) 
can be directly used to represent the kinetic energy of 
the rr-sector in Eq. ( p7|) in terms of c kx fermions: 



H r 



A L N 2 x /2 +Y,qb\ x b q 

q>0 



k 



k:c\ x c- kx :+A L P/& 



(53a) 
(53b) 



As a check, note that this equation also follows from the 
following observations: firstly, the equa tion of motion for 
the field tpx{x), express ed a s (44a) or (44b), is the same 
when calculated using (53a) or (53b), respectively, and 
therefore the latter two expressions can differ only by a 
constant; and secondly, this constant can be determined 
to be A^P/8, by requiring the free ground state energies 
for |0)s cxt given by the two expressions to be the same. 

Finally, in the subspace iS p h ys [of (p5|)] and hence also 
in S cx t, we can use (El) and (El) to express Af s S z and 



Af 2 in terms of c\c d - Thus, the EK-transformed 2CK 
Hamiltonian of Eqs. ([36]) and (|37| ) takes the form 



H' = H csf + H X + E G + const. , 



H, 



csf 



Y.J2i b lv b i 

c,s,f q>0 



H x =e d : c d c d : + ^ k : c\ x c kx : 

k 

+^AlT^2(cI x + c- kx )(c d - c d ) , 



(54) 
(55) 



(56) 



En = A; 



A4(l - Po) + (AA 2 + Aff + S 2 T - l/4)/2 



P/8 +e d {nf -1/2) + S T h e 



(57) 



The charge, spin and flavor degrees of freedom in H cs f 
evidently decouple completely. H x in (|56| ) has the form of 
a quadratic resonant level model whose "resonant level" 
has energy Ed and width T, where Ed = hi — h e is the en- 
ergy cost for an impurity spin-flip, and T = A 2 L /4a, which 
will be identified below as the Kondo temperature. 

Eq is the "free ground state energy" of the subspace 
S cx t in the presence of magnetic fields. Its SV^e term 
implies that the magnetic fields do not enter only in the 
combination hi — h e of s d , thus the role of the magnetic 
field h e applied to the conduction electrons is somewhat 
different from that of the local field hi. Note though, 
that for h e = 2nA^ (with n £ Z) the Sxh e term can 
formally be absorbed (up to a total energy shift) by in- 
troducing a "new total spin" S' T = St + 2n, since then 
A L S T /2 + S T h e = A L S T l /2 - 2n 2 A L . Now, since the 
construction of the complete finite-size spectrum involves 
enumerating all possible values of St, and since the gen- 
eralized gluing condition ( |75| ) to be derived below is in- 
variant under St ^ St + 2n, the finite-size spectrum for 
h e = 2uAl and a local field hi (so that e d = hi — 2nAi/) 
will be identical to that for h e = and a local field of 
hi — 2uAl (so that e d is unchanged). The physical origin 
of this "periodicity" is that as h e increases, at each value 
2uAl a "level crossing" occurs in which the free-electron 
ground state changes from, say, \M C , Af 8 , Aff, Af x )o to a 
new one differing from it only in the spin quantum num- 
ber, namely \Af c ,Af s — 2,Aff,Af x )o, by flipping the top- 
most spin-up electrons in both channels j ' = 1 and 2 to 
spin down. 

For general values h e ^ 2nA^ , there is no such symme- 
try (essentially since electron-hole symmetry in the spin 
sector is lost), and the corresponding finite-size spectrum 
differs from that at the periodicity .points in that some 
additional splitting of states occurs .Ell For simplicity we 
henceforth set h e = and consider only a local magnetic 
field, with s d = hi, but the more general case h e ^ can 
be treated completely analogously. 







V. FINITE-SIZE SPECTRUM OF 2CK MODEL 

EK studied the resonant level model H x of ([56|) in the 
continuum limit L — > oo. They mainly analyzed its im- 
purity properties, showing that these have NFL behavior 
because "half of the impurity" , namely Cd + c ., decou- 
ples. By keeping L finite, one can extend their analysis 
to include also the NFL behavior of electron properties. 
In this section, we illustrate this by diagonalizing H x at 
finite L and constructing its finite-size spectrum in terms 
of its exact eigenexcitations. 



A. Diagonalization of H x 

Since H cs f is trivial, we just have to diagonalize 
the resonant level part H x in the extended subspace 
Sext(Sr, A/" c , A//), which is straightforward in principle, 
since H x is quadratic. However, care has to be exercised, 
in particular regarding normal ordering: the change in 
ground state energy due to the interaction turns out to 
be of order — r, and the sub-leading (state-dependent) 
contributions of order relative to this energy have to 
be extracted carefully when constructing the finite-size 
spectrum. 

As first step, we define new fermionic excitations, illus- 
trated in Fig. @, whose energies are strictly non-negative, 



a- k = (c- kx + c f _- kx )/V2 1 
/% = -i (c kx - cl- kx )/V2 J 

for jfe = if P = 1, 
for Ed > 
for e d < ' 



a = ' : 

ad 



.t 

0.r 



,t 



for k > 0, (58a) 

(58b) 
(58c) 



and which have the virtue that the j3 k excitations decou- 
ple completely from the impurity: 



H x = ^2 ka l a k + E kf3 t^ + \ £ d\ a d a d 

k>0 k>0 

+ Y, V k( a l + a k)( a d~al) . (59) 



k>0 



Here the possible k- values are given by Eq. (45), and the 
hybridization amplitudes V% by 



V a = Vi 



fe#0 



2 = e l 



(60) 



Note that in (|58| ) we purposefully defined ajj and & such 
that the free reference ground state |0)s cxt , by (|50|), con- 
tains no ajj or excitations, i.e. 



a d |0) SeKt = a- k \0) s ^ t = /3 E |0) Sejrt = , 



(61) 



as illustrated in the middle entries of the first rows of 
Figs. |(a), |(b) and |(c). Note too that a^0) 5cxt is de- 
generate with |0)s ext if Ed = 0, as is ao|0)s oxt in the odd 
electron sector, P = 1. 

Since the Hamiltonian Eq. ( |5^ ) is quadratic, it can be 
brought into the diagonal form 



H * = £ * Pi?* 



k>0 



£>0 



5Er 



(62) 



describing a ground state energy shift SEq and non- 
negative-energy excitations relative to a reference state 
|0)s ext , an exact ground state of H' in <S ex t; defined by 



a e \0)i 



0k\O)i 







(63) 



This diagonalization can be accomplished by a Bogoli- 
ubov transformation of the form 



at = 



E 



ne{k,d} v=± 



B env {a\ + va n )/2 



(64) 



where the new operators a e are required to satisfy: 



e a 



[H x ,a 
{a\,a e >} = 5 £e > 



(65) 
(66) 

Eqs. ( p5| ) and (^) yield a closed system of equations 
for the eigenenergies e and the coefficients B env . In Ap- 
pendix |c| we solve them explicitly (by transforming to 
conveniently chosen Majorana fermions), with the follow- 
ing results. 

The excitation energies e are the non-negative roots of 
the transcendental equation 



e47rr 



cot7r(e/A L - P/2) 



and the ground state energy shift is 



5E C 



M 

2 



^ 2 

fe>0 



SI 

e>0 



For e > 0, the coefficients B env are given by 
B e d+ = g(s) \£d\ , B Ed - = g{e) e , 



B e - k + = £>(e) 



me 2 



t \ 2V k ek 

B ek- = ^)zr^ 



(67) 

(68) 

(69a) 
(69b) 



where the normalization factor g(e) is 

2A L r 



g(e) 



-,1/2 



±(e 2 - e 2 ) 2 + A L r(e 2 + e 2 d ) +47r 2 r 2 e 2 



(70) 



For e = 0, the coefficients Bo m m ust be considered sep- 
arately and are given in Appendix C 2 b. 

Eqs. (§§), @, @ and @ to®, together with 
the gluing conditions ( |75| ) discussed in the next subsec- 
tion, constitute a complete, analytic solution of the 2CK 
model along the EK line. 



10 



B. Evolution of Excitation Energies 



The eigenvalue equation (|67|) is a central ingredient of 
our analytical solution, since it yields the exact excita- 
tion energies e of H x , and also allows one to explicitly 
identify the various crossover scales of the problem. Let 
the label j = 0, 1, 2, . . . enumerate, in increasing order, 
the solutions Ej t p of ( |67| ) in a sector with parity P. Their 
smooth evolution as functions of T and \e^\ can readily 
be understood by a graphical analysis of Eq. (|67|), see 
Fig. ||, and is shown in Figs. ||(a) and ||(b) for P = and 
1, respectively. All but the lowest-lying (j = 0) solutions 
can be parameterized as 



£ 3,P 



At [j - | - f + S jtP ] 



.7 = 1,2,3,. 



(71a) 



where Sj,p £ [0,1] is the shift of e^p/Al from its 
r = Ed = value and is determined selfconsistently [from 
©] by 



1 1 1 

— I — arctan — 

2 7T 4-7T 



Jfr _ £j,p 

hp r 



(71b) 



with Th = e d /T. The lowest- lying modes are given by 



£o,o 







for e d = 0, 



(72a) 



A L { (-1/2 + Sq,q) e (0,1/2] fore^O, 
e ,i=0 foralir,e rf (72b) 



(see als o Ap pendix |C 2 b[ ) . 

Eq. ( |71b| ) shows very nicely that T and Th are cross- 
over scales: Firstly, in the absence of magnetic fields, 
i.e. for \Ed\ — \hi\ = T h = 0, the spectral regime be- 
low r is strongly perturbed [Sj.p ~ 1/2 for Sj t p <C T], 
whereas above T it is only weakly perturbed [6j t p ~ 
for Ej^p ^S> r]. This also follows directly from a graph- 
ical analysis of the eigenvalue equation (|67|). It is thus 
natural to identify the crossover scale T with the Kondo 
temperature, Tk — T. 

Secondly, in the presence of a local magnetic field, 
Th = hf/T > furnishes another crossover scale. When 
considering the T^-induced shifts in 5j t p relative to their 
values for Th = 0, se vera l cases can be distinguished [by 
direct inspection of Q71b| )] : 

(i) For Th <C Al, i.e. for \hi\ much smaller than a 
crossover field h c ~ \fTA~£, none of the T^-induced shifts 
are strong. 

(ii) For Th 3> Aj,, T, the crossover scale Th divides the 
spectrum into two parts: the TVinduced shifts are weak 
for all levels with e ^ Th, but strong for all those with 

(iii) For r > 1], > [a special case of (ii)] one 
can distinguish three physically different regimes: the 
spectrum is NFL-like (non-uniform level spacings) in the 
intermediate regime Th £ T, and Fermi-liquid 
like (with uniform level spacing) in the extreme regimes 
£ ^S> T and e -c Tfe. In the last of these regimes (right- 
most part of Fig. 0), the set of lowest-lying e's is in fact 



identical to that for the free case = 0, T = (leftmost 
part of Fig. |]), except that the free case has one more 
£ = mode, associated with the impurity level's two-fold 
degeneracy due to spin reversal symmetry for \hi\ = 
(which is broken if Th ^S> Al, implying \hi\ ^> A^). This 
shows quite beautifully that a magnetic field very liter- 
ally suppresses the effect of spin-flip scattering for the 
low-energy part of the spectrum; heuristically, this hap- 
pens, of course, since low-energy electrons do not have 
enough energy to overcome the Zeeman energy necessary 
for a spin flip in a magnetic field. 

Since at a finite temperature physical quantities are 
governed mostly by excitations of energy e ~ T, they 
will show NFL behavior for V 3> T^^Tp,^nd Fermi liq- 
uid behavior for T > V or T < T/JEoHeI, as sketched 
in Fig | 



C. Generalized Gluing Conditions 



A general eigenstate of H x in 6> cx t has the form 



Ma 



\E) ocJJ 4,114,10)^ 



(73) 



i=i j=i 



However, as emphasized earlier, of all such states only 
those in the physical subspace 5 p hy B must be retained, 
and all others discarded as being unphysical. (Recall 
that we had to extend 5 p h ys to <5> ox t to define the oper- 
ators ip x and Cd ) To identify which \E) are physical, 
we shall now derive a generalized gluing condition satis- 
fied by them, by noting that \E) can be physical only 
if the state \E) = limr^o \E), to which it reduces when 
r is adiabatically switched off, satisfies the free gluing 
conditions (|22|) . Key to the derivation is the fact that al- 
though the hybridization interaction H' ± of (|42]) does not 
conserve the number of a^ n excitations, it does conserve 
the parity of their number. 

To be explicit, let "Pg be the the parity of the number 
of excitations of \E) relative to |0)s oxt : 



(E\[ 



mod2|.E). (74) 



fc>0 



During the adiabatic switch-off of V, this quantity of 
course remains fixed, and hence equals P^(T — * 0). This 
in turn can be written as 



n=d,k>0 



k>Q 



= (E\ \ N x +a\a d \ moA2\E) 

= (E\ [(Af x -P)+tf a -Sr-% + 4 0) ] mod 2\E) . 

The first equation follows because the hybridization in- 
teraction preserves the parity of the excitation numbers 
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(in other words, since a\ is a linear combination of a^ n 
and a n ); the second follows because the c\ excitations 

kx 

counted by Sf x are linear combinations of ag, at, (3p. and 

/?| [this is illustrated by Fig. |], which depicts how states 
in the c^ x and aj., /% representations are related to each 

other]; and the third follows from (|5l]) for Af x and (58c), 
(E3) and (23) for ay. Imposing now the condition that 



\E) must be in the physical subspace tS P h ys and hence 
must satisfy the free gluing condition (|22b|), we obtain 



[N c + Nf 



S T -{P 
St-(P 



l)/2] mod 2 
l)/2] mod 2 



{e d > 0) , 
{e d < 0) . 

(75) 



This generalized gluing condition specifies which of all 
the possible states in S cxt are physical, i.e. are in iS p h ys ; 
it supplements the free gluing condition ( |22a| ), which stip- 
ulates that St ± 1/2 must be integer (half-integer) if M c 
and Aff are integer (half- integer). The unphysical states 
in (S ex t that do not satisfy (f75|) must be discarded when 
constructing the finite-size spectrum (and thus no such 
states occur in Table | below) . 



D. Ground State Energy Shift 

The form of Eq. ( |o8|) for the change in ground state 
energy SEq suggests that it can be interpreted as the 
dynamical binding energy of the impurity spin, which re- 
sults from the impurity-induced energy shifts of all the 
states in the fifed Fermi sea. (The factor 1/2 in ( |6§| ) re- 
flects the facttaH that only "half" of the x-pseudofermion 
field, namely ipx + ipl, couples to the impurity [in (ff2|)|.) 
For Sd — 0, the n umb er of levels strongly shifted by the 
interaction is [by (71b)] of order T/Al, and each of these 
gets shifted roughly by A^/2; we can thus estimate that 
|<5Pg| will be of order T, which supports the heuristic 
statement that "the impurity's binding energy is of or- 
der TV' . 

However, since the level shifts A^ <5 7 -.p a lso have a 
P-dependence of order ~ A|/r [from (71b)], the total 



ground state energy shift SEq will have a P-dependcncc 
too, of order ~ Al . We therefore write 



SEr 



SEZ + PSSZ, 



(76) 



where the first term is P-independent and hence gives 
only an overall energy shift. In contrast, SEq affects the 
finite-size spectrum since it shifts the odd electron states, 
P = 1, relative to even electron states, P = 0, and hence 
must be e valu ated with particular care. This is done in 
Appendix C4, where we find, for T/Ar, ^> 1, 



5E% = 



5E% 



-A L /8 (T h = 0), 
(T h » A L ), 

-2r[ln(P/4vrr) + l] , 
-2r[ln(D/| Ed |) + l] , 



{T h = 0), 

(T h » A L ,r). 



(77) 



(78) 



Here D ^> T, Th is a cutoff needed to regularize the sums 
in @. Note that for T h = 0, @ is consistent with the 
estimate for SEq above if we take D ~ 1/a and recall 
that T = \\fAa. For T h > T, the magnetic field \e d \ 
takes over as lower energy scale in the logarithm instead 
ofT. 



E. Construction of the Finite-Size Spectrum 



Now we are finally ready to construct the finite-size 
many-body excitation spectrum of the 2CK model. In 
doing so, we shall generally use calligraphic £'s to denote 
dimensionless energies measured in units of A^. Specifi- 
cally, we shall construct the dimensionless energies 



£{L) = [P(i)-P min (L)]/A L , 



(79) 



associated with the lowest few exact many-body eigen- 
states \E) of the full Hamiltonian H' of (p4[), measured 
relative to its ground state energy, P m in- For the sake 
of simplicity we only consider the case with periodic- 
ity index Pq = 1 [see (H)], for which the t/Jaj's have 
anti-periodic boundary conditions. In this case the free 
ground state in the electronic sector is unique, namely 
|0)o, which somewhat simplifies the counting of states. 
(Of course, one can use the same procedure for Pq = 0, 
with similar results.) 

The construction proceeds in three steps: we first 
evolve toward the EK line, second evolve along the EK 
line, and third turn on a local magnetic field. The re- 
sults are summarized in Fig. || and Tables | and |l[ The 
caption of Tabic | also summarizes the technical details 
of the construction. Here we just state the main ideas: 

(i) Phase-shifted Spectrum: — First we study the evo- 
lution of the spectrum toward the EK line for X z £ [0, 1] 
at Aj^ = Ed = 0. Since the impurity has no dynamics 
for Aj^ = 0, the spectrum is that of a free-electron Fermi 
liquid with a S z -dependent phase shift in the spin sector, 
given by iT'(A_i_ = 0) of (|3q); it evolves linearly with in- 
creasing A z , from £f rcc at A 2 = to £ p hasc at A z = 1, see 
Fig. |(a). 

(ii) Crossover Spectrum: — Next we study the spec- 
trum's further evolution along the EK line for T/Al e 
[0, oo) at X z = 1, Sd = 0. To this end we enumerate in Ta- 
ble | the lowest- lying physical eigenstates \E) of the full 
Hamiltonian H' in terms of the excitations at , fit and 

£ j,p 1 1 k 

b^y which diagonalize it, and follow the evolution with 
increasing T / Aj, of the excitation energies Ej t p (shown 
in Fig. |]) , and of the ground state energy shift SEq [see 
@]. This yields the crossover shown in Fig. ^(b) from 
the phase-shifted to the NFL fixed point spectrum, con- 
sisting of a set of universal, dimensionless energies defined 
by 
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£nfl = lim 



E(L; e d = 0, T) - E min (L; e d = 0, T) 



(80) 



Satisfyingly, the spectrum of £nfl energies found in 
Fig. ||(b) and Table [n] (degeneracies are given in brack- 
ets) coincides witiij— tte _anes obtained in NRG and 
CFT calculations M^mEl This constitutes a direct and 
straightforward analytical proof of the soundness of the 
latter approaches. In particular, it proves the so-called 
fusion hypothesis employed by Affleck|-aiid. Ludwig in 
their CFT calculation of this spectrum.E-TE3 As is well- 
known from CFTpil each of the fixed-point values <?nfl 
can be associated with the scaling dimension of one of the 
operators characterizing the fixed point. The occurrence 
of £nfl's that are not simply integers or half-integers is 
thus a very direct sign of NFL physics, since these corre- 
spond to non-fermionic operators. 

Our NFL spectrum demonstrates explicitly the well- 
known fact that the spin anisotropy is irrelevant at the 
NFL fixed point, since if we take the continuum limit 
Al — * at fixed T, the fixed point spectrum is evidently 
reached independently of the specific value of T. Said 
more formally: the symmetry of our anisotropic start- 
ing Hamiltonian with respect to transformations in the 
charge, spin and flavor sectors is U (l) c x U(l) s x SU (2)/, 
i.e. in the spin sector it is only invariant under spin rota- 
tions around the z axis; in contrast, Affleck and Ludwig 
derived the NFL fixed point spectrum by assuming it 
to have the complete U(l) c x SU(2) S x SU(2)f symme- 
try of the free model. The fact that the low-energy part 
(s <C Tk) of our NFL fixed point spectrum coincides with 
theirs beautifully illustrates how the broken symmetry of 
the original model is restored in the vicinity of the NFL 
fixed point, and thus proves another central assumption 
of the CFT solution of the 2CK model, in agreement with 
the NRG study of Pang and CoxEj 

The fact that the exact eigenenergies of H' interpolate 
smoothly between their values for Aj^ = and Ai ^ 
[Fig. ^(b)] may at first seem somewhat surprising, be- 
cause a common way of heuristically characterizing a 
NFL is that its quasiparticles are orthogonal to the bare 
ones of the corresponding free Fermi liquid. This is re- 
ferred to as the "breakdown of Landau's quasiparticle 
construction" , since in Landau's picture of a Fermi liquid, 
the dressed quasiparticles and the corresponding bare 
ones have finite overlap. Here, in fact, one can readily 
check that s cxt (0|a^.a^j |6)s e]rt is non-zero (where e(k) is 

the excitation energy that reduces to k as T/Al — > 0), 
implying that in the a-basis the system is a Fermi liquid 
(in accord with the fact that H x is quadratic, i.e. "sim- 
ple"). However, this does not contradict the fact that in 
the original Ckaj basis the system nevertheless behaves 
like a NFL, since the bosonization-refermionization rela- 
tion between states in the a and Ckaj bases is very highly 
non-linear. 

(iii) Crossover due to local magnetic field: — Finally, we 



turn on a local magnetic field, s d = hi =/= at fixed A z = 1 
and T/Al 3> 1, thus breaking spin reversal symmetry. 
The further evolution of the excitation energies Ej t p as 
functions of increasing Th/ Al, shown in Fig. |f|(b), yields 
the magnetic-field-induced crossover, shown in Fig. |^(c), 
from the NFL fixed point energies £nfl to a set of en- 
ergies £ p h corresponding to a phase-shifted Fermi liquid 
fixed point. For T^/A^ ^> 1, the impurity level evidently 
becomes empty for all low-lying states, (c d Cd) — 0, i.e. 
the impurity spin is frozen in the state S z = JJ- In- 
deed, the spectrum £ p h which one recovers is precisely 
the same phase-shifted spectrum as £ p hasc at the point 
A 2 = 1 and Aj^ = 0, apart from a degeneracy factor of 
two, due to the lack of spin reversal symmetry, compare 
Table || This shows very nicely that the magnetic field 
"erases" all traces of NFL physics for the lowest-lying 
part of the spectrum, because it effectively switches off 
spin-flip scattering. 



F. Finite-Size Behavior of Physical Quantities 

Let us now briefly discuss the finite-size, T = behav- 
ior of three physical quantities at the NFL fixed point, 
namely the entropy, susceptibility and the fluctuations in 

The entropy of the ground state at T = 0, e d = is evi- 
dently simply ln2 for any L, since the ground state is twp=. 
fold degenerate (see Fig. ^ . This should be contpasteoS 
with the celebrated result \ In 2 that one obtainsES when 
taking the limit L — > oo before T — > 0. The difference 
simply illustrates the triviality that the order of limits 
does not commute, since for finite L the system is always 
gapped. 

The susceptibility at T — due to a local field hi is de- 
fined by x = — a -)f ° ■ Since Eg = Eg + 5 Eg, we simply 



have to evaluate [by (|57|), (|68|)] the sum X = \ J2 6 fst 
For hi = 0, the summands can be determined by differ- 
entiating (pTj), giving 



x (h i = o) = J2l 



4ttFA l 



e>0 
1 



£ (A L 47rr + 7r[(47rr) 2 + e 2 }) 



(81) 



47r 2 r 



ln(47rr/A L ) 



(for f> Ai) . (82) 



The fact that x&i = 0) — > oo as i ^ oo is of course a 
characteristic sign of 2CK NFL physics: it illustrates the 
instability of thp-NFL phase with respect to a local sym- 
metry breakingH At finite temperatures T takes over the 
role of the infrared cutoff A^, so, that the susceptibility 
diverges logarithmically with TElla 

The fluctuations in J\f x can be quan tified by calcu- 
lating (N%) — {N x ) 2 - In Appendix C5 this is done at 
Ed = for the physical ground state of iS p h ys for both 
P = and 1. We find that (M x } = for arbitrary ratios 
of T/Al, showing that the ground state contains equal 
amounts of spin from both flavors j = 1,2, as expected 
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from the 2CK model's flavor symmetry. Furthermore, 
we find that (N x ) = P/4 for F/A — > 0, as expected intu- 
itively, since in this limit the considered ground states are 
linear combinations of states with M x — ±P/2 [compare 
Fig. ^|. In contrast, in the limit T/Al » 1 wc find 
that the fluctuations diverge logarithmically with sys- 
tem size, (N%) ~ ^j-lnTi. This illustrates very vividly 
how severely the dynamical impurity stirs up the original 
Fermi sea at the NFL fixed point. 



VI. CONNECTION TO DIFFERENT 
RENORMALIZATION GROUP METHODS 

In the literature several renormalization group (RG) 
methods have been applied to the multichannel Kondo 
model. In this section we relate these RG methods to our 
finite-size bosonization technique, by showing how the 
strategies employed by the former can be implemented, 
in an exact way, within the latter. 



A. High-energy cutoff scaling techniques 

The most common types of RGs are the ones used in 
particle physics and in the standard treatment of critical 
phenomena. In these RG procedures, one reduces a high- 
energy cutoff, say D, in order to gradually eliminate some 
high-energy degrees of freedom, arguing that they only 
slightly influence the low-energy physics of the system. 
The change in the cutoff must be compensated by rescal- 
ing the model's dimcnsionless coupling constants and 
masses in order to keep the physical properties (differ- 
ent inherent energy scales and dressed masses) invariant. 
These kinds of scaling procedures, which include Ander- 
son's poor man's scalingj£3 the multiplicative RGJIj and 
the Yuval- Anderson real time RG,Lj have been widely 
used in the continuum li^t_L-.pt. oo, to study the mul- 
tichannel Kondo model.y§Bli!fel 

In our case the high-energy cutoff D can be identified 
with the ultraviolet cutoff 1 /a introduced when defining 
the boson fields cf> a j [in (jljj)], i.e. D ~ 1/a. Thus, if a 
dimensionless coupling constant, say A, has the scaling 
equation 



dlnA 
d\nD 



dlnA 
din a 



= 7(A), 



(83) 



then in the weak-coupling regime A <C 1, its scaling di- 
mension is 7(0), and it is relevant, marginal or irrelevant 
for 7(0) < 0, =0 or > 0, respectively. 

Now, along thp_EK line one immediately obtains the 
scaling equationsEJ 



X 



din a 



A 2 = 1 



(84) 



The first, which follows from the requirement of the in- 
variance of the Kondo scale T = A^_/4a, shows that Aj^ is 



relevant and grows under bandwidth rescaling, wi th di- 
mension —1/2. As explained at the end of Section fVB, 
the sec ond eq uation follows from the fact, (proven in Sub- 
section VIC), that the leading irrelevant operator is ab- 



sent only at the EK line, where A 2 = 1. Eqs. (84) exactly 
coincide with the ones obtained with the Yuval- Anderson 
technique J23 




A± 



7T 



(85a) 
(85b) 



if in these the phase, shift 6 — A 2 7r/4 is replaced by 
7r/4 at th e EK line,E3 as discussed at the end of See- 
Remarkably, while these latter equations re- 
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tion 

suit from rather non-trivial calculations in terms of the 
original Hamiltonian, they can trivially be derived after 
the EK transformation, which in effect resums all appro- 
priate diagrams into a qua drat ic form. 

As discussed Subsection VB, in a finite local magnetic 



field Ed = hi, a further natural energy scale appears: 
Th = tij/T. For energies below this scale the magnetic 
field destroys the non- Fermi liquid behavior and a Fermi- 
liquid is recovered. By requiring the invariance of Th one 
immediately derives that, as long as the high-energy cut- 
off 1/a is much larger than the Kondo scale T, the field 
hi must be invariant under the RG transformation: 



dhi 
din a 



= (i/a>r). 



(86) 



However, once the cutoff is reduced sufficiently so that 
1/a -^i r, the role of T is taken over by 1/a, i.e. Th is now 
given by hf a, and therefore the scaling equation ^ 
must be replaced by 



dhi 1 , / , 

— — = --hi (i/a«r). 

dm a I 



(87) 



Since hi is dimensionfull, the scaling dimension of the lo- 
cal field can not be read off directly from this equation; 
instead one must consider the corresponding equation for 
the dimensionless field measured in units of the effective 
bandwidth, h = hid, namely 



(88) 



-— = -h (i/a«r). 

a m a I 

Close to the NFL fixed point the local field thus [by Q83 
has dimension —1/2 and is relevant: when measured in 
units of the effective bandwidth, it grows when the latter 
is decreased. (By a similar argument, the dimension of 
the local field in the regime l/a>T of ( |S6| ) is —1.) 

Equations (86) and (p7| ) are in complete agreement 
with those obtained by the Yuval-Anderson technique.EJ 
We remark at this point that perpendicular local mag- 
netic fields h x y Ii.e. perturbations of the form h x S x or 
hySy) are knowno to scale differently from hi — h z , and 
at the EK line their scaling dimension is known to be 
— 1/2 even in the region l/a>r. 
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B. Connection to Numerical Renormalization Group 

In this subsection we show that an analysis of our 
finite-size spectrum as function of L in fact represents 
an analytical version of Wilson's numerical renormaliza- 
tion groupEll (NRG), which can simply be viewed as a 
special type of finite-size scaling. 

In Wilson's procedure one divides the Fermi sea into 
energy shells using a logarithmic mesh characterized by 
a parameter A > 1, and then maps the model onto 
an equivalent one in which the impurity is coupled to 
the end of an infinite conducting chain, where the hop- 
ping between the sites n and n + 1 scales as A~™. The 
n'th site in this chain represents an "onion-skin" shell 
of conduction electrons, characterized by spatial exten- 
sion ~ A"/ 2 around the impurity site and energy ~ A - ™. 
The NRG transformation is then defined by considering 
truncated chains of length N with Hamiltonian Hn , and 
consists of (i) adding a new site to the end of the chain, 
i/jv — » ffjv+ij an d (h) rescaling the new Hamiltonian 
by A: fljv+i — > Ai/jv+i- Trivially, step (i) reduces the 
mean level spacing by a factor of 1/A, while step (ii) is 
needed to measure all energies in units of the new mean 
level spacing. This strategy is implemented by numeri- 
cally diagonalizing i?jv+i and retaining only the lowest 
few hundred levels. One finds that after a number of 
iterations the spectrum of Hn converges to a fixed, uni- 
versal set of energies, characteristic of some fixed point 
Hamiltonian.Ea For the 2CK model this spectrum has 
been shows to be identical to the one obtained by bound- 
ary CFT.E] 

The NRG procedure outlined above can easily be in- 
terpreted in terms of our finite-size calculations. Step (i) 
corresponds to increasing the system size, L — > AL, (i.e., 
reducing the level spacing, Al — > Aal = Al/A), while 
step (ii) is equivalent to measuring all energies in units 
of Aal. Combining both steps, an "analytical RG step" 
thus has the form: 



H x (L,T,e d ) H x {AL,T,e d ) H x (L,Ar,Ae d ) 



(89) 



where the last equality follows identically from Eq. (56) 



For Ed = this means that increasing the system size at 
fixed r is equivalent to increasing L at fixed L, empha- 
sizing once more that in this case the spectrum depends 
only on T/Al- Therefore the "spectral flow" as function 
ofr/A L in Fig. | can be viewed as the analytical version 
of an NRG spectrum as a function of iteration number. 



C. Finite-Size Scaling 

It is also straightforward to implement Wilson's pre- 
scription for extracting the exact scaling exponent of a 
perturbation around the fixed point, say 5XO, from its 
effect on the finite-size spectrum: In general, it causes the 
dimensionless energy £~{L) [of J79|)] (calculated at a finite, 



non-zero Al <C L) to differ from its universal fixed point 
value £nfl [of (|80|)] by an amount S£(L), whose leading 
asymptotic behavior for L — * oo is 



5£{L) = £{L) - f nfl ~ (6X/L^) n 



(90) 



where n > 1 is some integer. By definition, 7 is the scal- 
ing dimension of the operator O, which is relevant for 
7 < 0, marginal for 7 = and irrelevant for 7 > 0. Thus 
deviations from the universal spectrum are characteristic 
of the operator content of the fixed point. 

We first consider the situation on the EK line (i.e. for 
X z = 1), and close to the NFL fixed point, where Al/F 
and Th/ Al are are both 1 (at NFL fixed point they are 
both 0). The leading deviations e^pj Al — (e^p/A^nfl 
of the dimensionless single-particle excitation eigenener- 
gies from their NFL fixed point values are then given 
[from @] by 



1 



S j,p - (<5j,p)nfl - -r-n 



Al(j-P/2) 



4^ 2 [A L (j-P/2) 



(91) 



(for j > 1; for j = 0, Eq. ( |C18| ) yields the same con- 
clusions) . The leading dependence on the local magnetic 
field via T h /A L is evidently [h l L l/2 ] 2 , which grows as 
L — > 00. This shows that a local magnetic field has di- 
mension — — 1/2 and is relevant: for an arbitrarily 
small hi, there exists a system size L above which the 
lowest part of the spectrum and the ground state prop- 
erties of the model are drastically affect ed, n amely when 
Al < Th) in the language of section VB, this occurs 
as soon as the crossover field below which the spectrum 
is unaffected, namely h c = VLAl, becomes smaller than 
\hi\ as L is increased. The dimension 7^ = — 1/2 also fol- 
lows from the L-dependence of h c , and agrees with the 
conclusions of our bandwidth rescaling arguments [see 

In the absence of magnetic fields, the leading term in 
j9l| ) vanishes with increasing L as (TL)^ 1 , implying that 
the least irrelevant irrelevant operator on the EK line 
has dimension Jek = 1- Thus, we conclude that the 
leading irrelevant operators with diiaension 7=1/2 that 
were found in the CFT treatmentE3 are absent on the 
EK line, in agreement with Refs. |^^p|. 

Now let us move away from the EK line by taking 
A 2 = 1 + 5X Z , and do perturbation theory in SX Z , i.e. in 
SH' Z of (|3^). Then the operators with dimension 7=1/2 
just mentioned immediately show up: As shown in de- 
tail in Appendix |C q , we find that the "zero mode" term 
SX Z AlNsSz of Q38|) (which does not occur in the con- 
tinuum limit considered in Ref. ^), affects the spectrum 
already in first order in 6X Z : in the absence of magnetic 
fields, the first excited states (with £nfl = 1/8) are 
shifted relative to the doubly degenerate ground states 
(with £nfl = 0) by an amount 



S£(L) ~ -|*A,(1 + 4^ 2 r/ALr 1/2 ~ L-^ 



(92) 
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This implies that the leading irrelevant operator away 
from the EK line has dimension 1/2. 

In the presence of a local magnetic field Ed = hi, one 
finds in the continuum limit A^ <C T, hi that the ground 
state degeneracy is split by an amount 



SX Z \hi 



In- 



S£(L) = 



2ir 2 T 4vrr 



sx z 



4r 



(A L < hi < r) , 



(A L < r < 



(93) 



This shows that the magnetic-field behavior along the 
EK line is not completely generic, since it misses this 
part of the /^-dependence of the magnetic-field-induced 
crossover. Note that the | hi \ /T In |/ii|/r behavior that 
occurs for a local magnetic field of intermediate strength 
is consistent with the conclusions of the NRG studies of 
Rcf . ^ for the /^-dependence of a certain phase shift that 
can be used to characterize the NRG spectra. 

Finally, we would like to comment here on the iden- 
tification of the Kondo scale Tk- In Section VB we 



showed that the crossover scale below which the finite- 
size spectrum takes its fixed-point form (at hi = 0) 
was T, and hence concluded that Tk — T. This dif- 
fers from the suggestion of Sengupta and Georgealj that 
the Kondo scale in the anisotropic 2CK model close to 
the EK line is not T but rather T/(SX Z ) 2 . This scale 
emerged naturally in their calculation of the total suscep- 
tibility enhancement due to the impurity, which yielded 
Ximp ~ (<SA 2 ) 2 /Tln(r/T) (at h, = 0). However, the fac- 
tor (6X Z ) 2 only expresses the fact that the amplitudes of 
the leading irrelevant operators vanish on the EK line, so 
that the characteristic logarithmic features appear only 
in second order in SX Z . The fact that the scale above 
which these logarithmic features vanish is T ~ T, not 
T ~ T/(SX Z ) 2 , supports our above conclusion that it is 
rather T that should be identified as the Kondo scale. 



VII. DISCUSSION AND CONCLUSIONS 

The main general conclusion of our work is that con- 
structive finite-size bosonization is an unexpectedly pow- 
erful tool for investigating quantum impurity problems. 
Firstly, for the 2CK model, it enables one to analyti- 
cally calculate by elementary means the crossover along 
the EK-line of the finite-size spectrum (and the corre- 
sponding eigenstates) between the free Fermi liquid and 
the NFL fixed point. This crossover had hitherto been 
tractable only with the numerical renormalization group, 
and has been beyond the reach of all analytical ap- 
proaches used to study this model. 

Secondly, finite-size bosonization can deal without 
much additional effort with symmetry-breaking perturba- 
tions, such a&-a finite magnetic field (or channel symme- 
try breakingo, which was not discussed here, but can 
be included by a straightforward extension of our meth- 
ods). Indeed, it is to be expected that the methods de- 



veloped here can fruitfully be applied to a number of 
related quantum impurity problems. For example, an 
adaption of our finite-size refcrmionization approach was 
very recently used to rigorously resolve a recent contro- 
versy regarding the tunneling density p£ states at the site 
of an impurity in a Luttinger liquid.Cj Other potential 
applications would be to the generalized Kondo-niodels 
studied by Yeea or by Moustakas and Fisher&El, or by 
Kotliar and SiEJ. 

Thirdly, finite-size bosonization allows one to mimic 
in an exact way the strategy of standard RG approaches 
such as poor man's bandwidth rescaling and finite-size 
scaling; thus it should be useful also as a pedagogical 
tool for teaching and analytically illustrating standard 
RG ideas. 

Crucial to the success of our method is that we 
do not use field-theoretical bosonization, in which the 
bosonization relation just has the status of a formal cor- 
respondence, but Haldanc's constructive formulation of 
bosonization, in which all operators and fields needed 
are constructed explicitly in terms of the initial set of 
electron operators Ckaj in terms of which the model is 
defined. This has the great advantage that Emery and 
Kivelson's bosonization-refermionization mapping of the 
model onto a quadratic resonant level model can be im- 
plemented rigorously, not merely as a formal correspon- 
dence, but as a set of operator identities in Fock space. 
To achieve this, however, the Klein factors have to be 
treated with due care. Our main technical innovation was 
to demonstrate how refermionization can be performed 
at the same level of rigor as bosonization, by extending 
the Fock space to include unphysical states, and identi- 
fying and discarding these at the end of the calculation 
using a generalized gluing condition. However, we wish 
to emphasize that our method is truly elementary: in 
principle it requires nothing more than a knowledge of 
standard second quantization, since that is all onejieeds 
to derive the constructive bosonization formalisms 

Our rigorous implementation of EK's mapping of the 
2CK model onto the resonant level model ensures that 
the latter is not merely an "effective model that captures 
the essential physics" , but truly identical to the original 
model one starts out with (after the standard reinterpre- 
tation of the meaning of the coupling constants discussed 
in Appendix [A], to compensate for differences that can 
arise relative to other methods due to the use of differ- 
ent regularization schemes). The fact that we are thus 
able to solve the original model exactly, constitutes a sig- 
nificant advance relative to a number of alternative ap- 
proaches that have tried to analytically access the NFL 
fixed point, but end up doing so using either an "effec- 
tively equivalent" model or some assumptions that are 
not proven from first principles. Let us briefly mention 
three of these: r~K~~i 

(i) Coleman and coworkersE§E3 have proposed a 
"pedestrian solution" of the 2CK model, in which it is ar- 
gued that many of its properties can be calculated using 
a so-called "compactified model" involving only a single 
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channel of spinful conduction electrons. This model was 
argued to represent that part of the 2CK model that 
is left over when one "factorizes out" the charge and 
flavor degrees of freedom. Indeed, using field-theoretic 
bosonization, Schofield showed that there is a formal cor- 
respondence between the compactified model and our H±_ 
of Eq. (|2^) (which involves only ip s and <p x ), and that it 
yields the same results as the 2CK model for the impu- 
rity contribution to thermodynamical properties. In this 
sense, the compactified model can be viewed as an effec- 
tive model for calculating impurity properties. However, 
as first emphasized by Ye,t3 it is not equivalent to the 
original 2CK model, since Schofield's arguments ignored 
the fact that there are gluing conditions such as (g2|) be- 
tween the c,f sectors and the s,x sectors. As long as these 
are ignored, the compactified model can not be used to 
calculate conduction electron properties, since that re- 
quires adding back the contributions from the charge and 
flavor channels- This was attempted by Zhang, Hew- 
son and Bulla,E3 who calculated the fixed-point spectrum 
of the compactified model using the NRG, from which 
they tried to reconstruct that of the 2CK. However, their 
construction has an ad hoc character (it requires knowl- 
edge of the answer) and succeeded only partially (it did 
not correctly reproduce all degeneracies of the 2CK fixed 
point spectrum). 

Our constructive bosonization approach allowed us to 
clarify this issue completely: it makes precise in what 
sense the c and / sectors can be "factorized out", rig- 
orously yields an appropriate model for the remaining s 
and x sectors, emphasizes the gluing conditions between 
the c,f and s,x sectors, and shows how they can be used 
at the NFL fixed point to combine the contributions from 
all four sectors to obtain the NFL fixed point spectrum. 

(ii) Affleck and Ludwig's path-breaking conformal field 
theory solution of the model in the NFL regime, though 
very elegant and highly successful, rests on two assump- 
tions that can not be proven from within their theory, 
but have to be confirmed by comparison with other ex- 
act methods: firstly, they assume that the NFL fixed 
point has the same U(l) c x U(2) s x SU(2)f symmetry 
as the free model; and secondly, they need to use a cer- 
tain conformal fusion hypothesis to obtain the operator 



content of the NFL fixed point. In Section [v|, we ana- 
lytically proved these assumptions in a simple, natural, 
direct manner. (The second assumption has also been 
proven by Sengupta and George using bosonization, but 
their proof requires extensive knowledge of CFT, and the 
assumption that the fixed point is invariant under mod- 
ular transformations). 

(iii) Maldacena and Ludwig@ have used CFT to show 
that Affleck and Ludwig's CFT solution can be reformu- 
lated in terms of free boson fields ip y (^satisfying certain 
asymptotic boundary conditions. Yal3 reproduced this 
result using field-theoretic bosonization (in the contin- 
uum limit) and scaling arguments, which, however, have 
the standard weakness of scaling arguments in strong- 
coupling problems: they give the initial flow of the cou- 
pling constants, but become non-rigorous once the weak- 
coupling regime is left. 

We have shown in_Ref. [l5] (and will elaborate this in 
a future publication^) that these results can be repro- 
duced with great ease and much more rigor by simply 
taking the continuum limit L — > oo of our above finite- 
size calculation. In fact, this allows us to fully reproduce 
all Affleck and Ludwig's results for electronic correlations 
functions. 

In summary, we believe that our finite-size bosoniza- 
tion approach is the first straightforward analytical cal- 
culation which, starting from first principles and without 
any assumptions, yields the crossover of the finite-size 
spectrum of the 2CK model from the free to the NFL 
fixed point and allows a detailed finite-size scaling anal- 
ysis of the fixed points. 

Acknowledgments: — We are deeply indebted to 
Michele Fabrizio for carefully reading the manuscript 
and for invaluable help in the early stages of this work. 
We are also grateful to Andreas Ludwig and Abraham 
Schiller for illuminating discussions and generous help, 
Gabriel Kotliar and Anirvan Sengupta for useful sugges- 
tions, and Volker Meden for a discussion on cutoff-related 
matters. This research has been supported by SFB195 of 
the Deutsche Forschungsgesellschaft and by the Hungar- 
ian Grants OTKA T026327 and OTKA F016604. G. Z. 
has been supported by the Magyary Zoltan Scholarship. 



APPENDIX A: CUTOFF-RELATED MATTERS 



1. Extending the Bandwidth to Infinity 



In this appendix we discuss in some detail various mat=. 
ters related to the choice of an ultraviolet cutoff scheme ,E3 
since this is a rather subtle matter, which can be elu- 
cidated more explicitly when using constructive rather 
than the more usual field-theoretic bosonization formal- 



ism. In Section A 1 we explain the need to reintroduce 
an ultraviolet cutoff after remov ing the fermionic band- 
width cutoff D, and in Section A 2 discuss how this is 



accomplished by the bosonization cutoff scheme used in 
the main text. 



Onphysical grounds, the momentum sums in Eqs. (Q) 
and (js ) for and i?; n t must be cut off at some large mo- 
menta (\k\ < D ~ pp) , to account for the finite width of 
the fermion conduction band and the fact that a realistic 
impurity potential always has non-zero range. However, 
the bosonization procedure used by us requires a single- 
particle Hilbert space with an unbounded fermion mo- 
mentum spectrunuo To achieve this we have removed, 
following Haldanecj, the implicit cutoff D in Eqs. (jl]) 
and (||), and used k 6 [— oo,cx}] instead. By doing so, we 
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extended the Hilbert space of single-electron states to in- 
clude unphysical "positron" states with arbitrarily large 
negative momenta, but this should not change the low- 
energy physics of the system, since by construction they 
require very high energies (> Ep) for their excitation. 

In the resulting single-particle Hilbert space with an 
unbounded fcrmion momentum spectrum, the fermion 
fields defined in ([}]) and the corresponding densities do 
not depend on any cutoff parameter. In this Appendix 
we shall denote this fact by a superscript (°> : 



r E 



(0), 



(Al) 



(A2) 



In this notation, the position-space representation of the 
Hamiltonian, given here to facilitate comparison with 
field-theoretic treatments of the 2CK model, reads: 



L/2 



2n 



H h = htS z + 



E V*,. =^(0)^^(0) 
, f L/2 dxa 



I -L/2 



2tt 2 



(A3) 
(A4) 

(A5) 



Now, some physical quantities, such as the phase shift 
5 of the outgoing relative to the incoming fields, 



tf£?(0-) S e«V$(0+), 



(A6) 



depend explicitly on an ultraviolet cutoff and in fact 
would be ill-defined without any. Therefore, the deci- 
sion to use an infinite fermion band must always be ac- 
companied by a reintroduction, in some other fashion, 
of an ultraviolet cutoff. Moreo v e r ,., th e precise way in 
which this is done is well-knowr£3EjO to strongly influ- 
ence the meaning of the coupling constants in i?i n t. As 
an example, we consider the case of no spin-flip scatter- 
ing (A_l = 0), and reintroduce an ultraviolet cutoff by 
replacing \ z in Eq. (B) for /£ n t eitherEJ by the separa- 
ble form A z ie - d fc l + l fc f) a / 2 , ora by the nonseparable form 
A Z 2e _ ' fc ~' c ' a / 2 . Choice 1 restricts both momenta of a 
fermion scattering process separately to a band of width 
1 /a, choice 2 only the momentum difference. They imply 
two different versions for H z in position space, namely 



2 

aj 



L/2 



dx 



' —L/2 



L/2 



L/2 



dx' S a/2 (x)S a/2 (x') 



(A7) 



H z2 = Xz2S. 



2 

aj 



L/2 



dxS a/2 (x):^(x)^(x):, 



-L/2 



(A8) 



as follows by noting that e ' fc ' a is the Fourier-transform 
of the smeared S a (x) function of ([lfjj): 



(A9) 



1 J2 e ~ikx e -\k\a = 5 a (x)+0(l/L 2 ). 



n k ez 



Eq. (|A7|) shows that choice 1 separately smears out both 



^a]\ x ) an d ^aj ( x ) over a range o,, i.e. corresponds to 



,(0) 



a zero-range potential in a finite band, whereas by (AS) 



choice 2 smears out the density Pa]{x), i.e. corresponds 
to a finite-range potential in an infinite band. Conse- 
quently, the equations of motions differ: Choice 1 yields 



i(d t -d x )i> ( °)(t,x) 



L/2 



dx'S a/2 (x')^(t,x') (A10) 



-aA-|.S". S a / 2 {x) . „ 

-L/2 



1 



(All) 



(we took a — » in the second line). This is solved byEl 

iP^(t,x) = e - lk ^ [9(x) + 0{-x)e 2iSl ] , (A12) 
§1 = — arctan(7raA 2 iS , 2/2) , (A13) 

where 0(x) is a sharp step function, and the phase 
shift agrees with that found in the Bethe Ansatz Kondo 
literature, pQr in non-ID treatments that use a finite 
bandwidtf£j. In contrast, the equation of motion for 
choice 2, 

i(d t - flt)^(t, a:) = na\ z2 S z S a/2 (x)^ (t, x) , (A14) 
has the solutionS 

l/>Vj(t,x) = e -«(*+*) e <aA. a S.ar C tan(!te/a) ( ( Al5 ) 

so that (for a < 1 =*= | ) the phase shift of (|A|) is 

S 2 = -na\ z2 S z /2 . (A16) 

Evidently, regularization schemes 1 and 2 yield differ- 
ent relations between coupling constant and phase shift. 
Since the latter, being a physical quantity, must have 
the same value in both schemes, Si = S 2 , we thus con- 
clude that the coupling constants must be related by 
A Z 2 = | arctan(7rA z i/4). 

Finally, note that after the removal of the fermion 
band cutoff D, even the free theory (ifi„t = 0) re- 
quires the reintroduction of an ultraviolet cutoff: the free 
imaginary-time-ordered zero-temperature correlator, 



(T^ {t,x)& (0,0)) = 



IX 



(A17) 



has a divergence at t = x = 0, which is often regularized 
by the replacement r — > r + sgn(r)a, where a ~ 1/pF 
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(though we reuse the notation a here, this cutoff param- 
eter in general need not be the same as that used in H z 
above). Alternatively, one can introduce a bandwidth 
cutoff into the definition of the fermion field itself, e.g. 

replace (x) by 



W s VTE e_fce " l * l ' /a ^ 



(A18) 



L/2 



dx'S a/2 (x - x')^(x') , (A19) 



L/2 



i.e. by a smeared version of ip^j (x), which results in 



t + ix + sgn(r)et 



(A20) 



2. The Bosonization Cutoff Scheme 

The bosonization cutoff scheme used in the main text 
constitutes yet another way, alternative to those just dis- 
cussed, of reintroducing an ultraviolet cutoff after re- 
moving the fermion bandwidth cutoff D: one bosonizes 
the theory completely and introduces an ultraviolet cut- 
off e~ 9 ' a '/ 2 in the boson fields of Eq. ([IJ), whose a- 
dependence we shall indicate in this appendix by an ex- 
plicit superscript ^ a ^: 



q>0 



51 (< 



Jqaj 



(A21) 



Evidently, 1/a can be viewed as "effective boson band- 
width" for the particle-hole excitations occuring in cf>( a \ 
In this notation, the bosonization relations ( p"8| ) and ( p~9| ) 
used in the main text read 

^ a a hx) = F aj A L /2 e- l ^- p °/ 2 ^ L ie-^W; (A22) 



a] 
» 



1 



(A23) 



[un- normal-ordering the exponential in (A22), which is 
possible only for a ^ 0, yields a factor (aA^) -1 / 2 (see 
Eq. (42) of Ref. ||) and thus reproduces ©]. The su- 



perscripts emphasize that for a ^ 0, rjj [a> and p (a> are 



no£ identically equal to the original ?/^ ' ) and p^ of ( Al) 
of (|A2|), which do not depend on a. Instead, the rigor- 
ously exact bosonizat ion id entiti es fo r ip^ and are 
the a = versions of (|A22|) and (|A23|) (cf. Ref. [5[ foot- 
note 7). Correspondingly, the exact bosonized position- 
space version of Hq depends on tp(°> (not 4>^): 



H = J2 A LN aj (N aj + l-P )/2 

aj 



By taking a ^ in ( |A22[ ) and ( |A23[ ) (but not in # ), as 
we do in the main text, we thus effectively make the re- 
placement ip^ — > ip( a >, i.e. we redefine the fermion fields 
to be explicitly cutoff dependent and thereby modify the 
ultraviolet behavior of the theory. Although this redefi- 
nition is not identical to the replacement ip^ — ► ip^ of 
(A18) since ip^ ^ tp^ a \ it similarly regularizes Mae cor- 
relation function [TijA a ^^ a ^) , which turns outo to be 
given by ( A2C| ) too. More over , it smears out the density 
by a, since [by ( |A2l| ) and @] 



p^(x) = A L J2e- lk - k ' la/2 :c{ aj c k , a -. 



(A25) 



L/2 
-L/2 



dx' 6 a/2 (x-x'): < )f (x')< } (x') : , (A26) 



and i/ius also regularizes H z , which depends on the spin 
density (cf. (|2£ ) for i? z in the main text). In fact, 
comparison of (A26) with H Z 2 of (A8) shows that our 
bosonization cutoff scheme regularizes H z precisely ac- 
cording to the choice 2 discussed above. 

It is therefore not surprising that the phase shift for 

[Vg via the EK- 



ip£j found at 



the end of Section 



o tA*S*¥>l a) (0) 



transformation U 
X z = 1, agrees with ( Alt ) for choice 2. 



namely |(5| = 7r/4 for 
(Note that U 

would be undefined for a ^ 0, since its exponential, in 
order to be unitary, must be non-normal-ordered.) How- 
ever, if one examines the phase shift more closely than for 
|x| ^> a, one discovers that the phase factor arctan(2a;/a) 
in (A15), obtained by solving the equation of motion 



for ip^J ' , differs from the arctan(x/a) in (|35|), obtained 
by EK-transforming ip^j • This simply illustrates that 
ip(°) ^ ,/,(<»). Indeed, if one EK-transfor ms ip ^ instead 
of V> (a) , one recovers the arctan(2x/a) of ( |A15| ), either by 
using i)£](x) cx e -i°«P?H*)/2 and 

U(pf ] {x)U^ = ip s (x) - 2X Z S Z arctan(2ir/a) , (A27) 

or by using the fermionic definition ( |A1[ ) for ip^) to find 

[<¥><»> (0), (x)] = ai arctan(2x/a) (x), (A28) 

together with the fact that if [A, B] = cB with [c, A] = 
[c, B] = 0, then e A Be~ A = Be. . 

This example illustrates the subtle difference between 
our bosonization cutoff scheme, which replaces ip^ by 
-0( a ) and thereby modifies H z — > H z2 , and the regulariza- 
tion scheme of choice 2, which modifies only H z — > i? Z 2 
but does not change ip(°\ As far as H z is concerned, 
both schemes can be used with equal merit, but once 
one has chosen one of them, one must use it consistently 
throughout. 

For the treatment of H±, however, our bosonization 
cutoff scheme is distinctly more convenient. To see this, 
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note that the EK-transformed version of H±_ differs, de- 
pending on whether H± is expressed through (as in 
the main text), or through In the former case the 

factor multiplying (X±/2)F^F- a j in (37) is 



APPENDIX B: CONSTRUCTING A BASIS FOR 
THE EXTENDED FOCK SPACE 



U 



A L :e iQ [^ a, (o)+M a, (o)] . g_ 

^y /2 :e -M»'(o) :5a) 



(A29) 



(for X z = 1), the latter case instead yields 



U 



A L : e M [^ 0) (°)+M 0) (o)]. S- c 

1/2 



LP 



= 2 " :e ia[ *'- 0)(:o)_a ' p - a5(0)+3 ' v '- )(0)1 : S- a . (A30) 

Whereas in the former case the (/^-dependence conve- 
niently drops out for arbitrary a, in the latter it in- 
conveniently drops out only in the limit a — > 0, in 
which the prefactor a -1 / 2 diverges (and moreover the 
fermion correlation functions are ultraviolet divergent). 
Note also that the extra prefactor of 2 in the latter 
case (which stems from normal-ordering the product of 



In Section IV C 1 we transformed from an old to a new 
set of quantum numbers, N — > A/", and embedded the 
physical Fock space .Fphys in the extended Fock space 
F cxt [ all \Jy ) e . Fphys satisfy both the free gluing condi- 
tions (22a) and ( |22b| ), but only ( [22a| ) is satisfied by all 
\M) G F CK t\- In this Appendix we show explicitly how 
such an embedding can be accomplished, by construct- 
ing a basis of A -particle ground states {\N )o} that spans 
^"ext, in terms of ordered products of new Klein factors 
Tly acting on two reference states. 

To begin, we fix the relative phases of the set { |iV) o } 
of N particle ground states that span F p h ys , by defining 



\N)o 



r 11 r .11 - r T2 .12 l u /0 



;,tJV T2 ; 

|2 



JN l2 | 
f2 



(Bl) 



States with an even or odd total number of particles have 
P = 2A r c mod2 = or 1, respectively. Clearly, all even or 
odd states can be generated, respectively, from the even 



or odd reference states |0)o or 



= pt 



F^\0) o , defined as 



US- a W = 



'MS-a and 



(°) : instead of 



: e 



ia(p ( 3 a} (0) . 



), implies that the coupling constant must be Klein factors F^ or F 



I f ) = \fi = f, f , f , f ) = \N = P, 0, 0, 0) , (B2) 
by the application of a product of an even number of old 



reinterpreted such that 2X± of the latter case corresponds 
to the Aj_ used in the former. This illustrates once more 
how sensitively the meaning of the couplings depends on 
the choice of regularization scheme. 



3. Point-Splitting vs. Normal-Ordering 

In the liter ature the position-space versions of Ho, 
Eqs. QA3D or (A24), are used more frequently than the 
momentum-space versions of Eqs. (Q) and (Q), perhaps 
because the former may seem more concise. The product 
of two fields at the same point is then regularized using 
the point-splitting prescription 



M/ ... . a j. By using Eqs. ( ffO| ) and related 
bilinear relations, this product can be transcribed into a 
product of an even number of new Klein factors T^ y or 

T y . The resulting state evidently is an eigenstate of J\f y , 
and since Eqs. ( [io| ) by construction respect Eq. (21), its 
eigenvalues N are related to N by Eq. ( pl|) . Therefore 



Wo 



_ p i*(Af) 



|A0o, 



(B3) 



where the Af particle ground state \N)o is defined to be 



|A7) ^t^t^t^t^|| )o 



(B4) 



lim 0\{x — ixq)02(x) 



(0\dx(x - ix )d 2 (x)\0) 



which in most cases is equal to the normal-ordered prod- 
uct : Oi{x)C>2(x) :, evaluated by normal-ordering the 
Ckaj's in the Fourier expansions of these operators (see 
e.g. Appendix G of Ref. E4T). However, when using point- 
split operators, great care is required if terms of order 
Ai are to be treated correctly. Since in practice they 
are more easily dealt with by using normal-ordering in 
the momentum-space representation than point-splitting 
in the position-space representation, we use the former 
throughout this paper. 



for A7 e (Z + P/2) 4 , the integers N y are defined by 
N y = My — P/2, and N y is, by construction, an even 

number. The phase factor e 1 *^- 1 = ±1 can be deter- 
mined, if necessary, by explicitly rearranging the above- 
mentioned even product of new Klein factors into the 
standard order of ( B4 ) . It ensures that the action of cor- 
responding pairs of old or new Klein factors on the l.h.s. 
or r.h.s. of Eq. (B3), respectively, produces the same re- 
sult. 

Evidently, the set {|A/)o} of all states with Y^y^y = 
even constitutes a basis for the physical Fock space ^"phys, 
just as {|-/V}o} does. The unphysical part of the extended 
Fock space J- cyL t can now be formally constructed by using 
the definition (B4) also for integers with J^Ny = odd. 
Note once again that t hese new states violate the second 
free gluing condition (22b) and are purely mathemati- 
cal constructs outside the original Fock space. Then the 
total extended Fock space can be formally written as 
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^^^^({|0)}®{||)}) 



where {|0)} and {|^)} denote the set of all states that 
can be generated from the reference states by the action 
of bosonic excitations b' qy . 

Within J-cxt , which is the natural extension of the orig- 
inal physical Fock space J-ph ys , the action of arbitrary 
(even and odd) products of new Klein factors evidently 
is trivially defined, and they leave the subspaces gener- 
ated by |0)o and |i)o separately invariant. Note though, 

that it is impossible to reach |0)o from |i)o or vice versa 
using new Klein factors, simply because these change N y 
by ±1, whereas the eigenvalues in the two subspaces 
"differ by i". However, they are of course connected by 

the original Klein factors, e.g. ||)o = ivJOjo- This shows 
again that there is no way to express an individual old 
Klein factor in terms of the new ones, or vice versa. 



APPENDIX C: EXPLICIT DIAGONALIZATION 
OF Hx FOR 2CK MODEL 

In this Appendix we diagonalize the Hamiltonian H x 
of ( |59| ) in explicit detail. We also calculate the ground 
state energy shift SEq, the Af x fluctuations (Af x ), and do 
perturbation theory about the EK line. 



1. Introducing Majorana Fermions 

Our aim is to find the unitary transformation that 
brings the Hamiltonian H x of ( |59| ) into the diagonal form 
(|62|), and to determine the discrete set of eigenenergies e. 
This transformation will map the original set of orthonor- 
mal operators occurring in (|59|), {a n ; n — d, k > 0} (with 

{anjttj,/} = 5 nn >), onto the new orthonormal set occur- 

=r,t 



ring in Eqs. ( p2| ) to (|66|), {a e ;e > 0} (with {5t e ,a\,} — 
<W); however, the transformation does not involve the 



,%'s in (59) at all, since they are completely decoupled 



and "just go along for the ride" below. 
Since the hybridization term in Eq. i 



the combinations 



(4 



aj,) and (ay 



9|) only contains 
"half of the 



impurity", (ad + ct^). is completely decoupled from con- 
duction electrons if £d — 0. (EK were the first to em- 
phasize that this causes the model's NFL behavior.) To 
exploit this fact, it is convenient to transform the two 
sets of fermions, {a n } and {<5 £ }, to two sets of Majorana 
fermions {^nv} and {lev} (f = ±), respectively: 



7rH 
Tn- 

''%+ ' 



1 1 

—i i 



1 

-( 1 1 



{n = d,k>0), (CI) 
(e > 0). (C2) 



By construction they are real (j nL , = 7^„, 7 £ „ — 7^) and 
obey the anticommutation relations 



(C3) 

(C4) 



When rewritten in terms of these Majorana operators, 
the original form (E9J) for H x becomes 



H , 



J2 Hm+Tk- + 1/2) + M(*7d+7d- + 1/2) 

fc>0 

k>Q k>0 



(C5) 



[with k and Vj, given by Eqs. (|45| ) and (pCj)], and the 
sought-after diagonalized form ( |62D becomes 

H x = J^s (* %+%- + 3) + SE G + X] ■ ( C6 ) 



e>0 



fc>0 



To find the transformation that brings (C5) into the form 
( C6 ) , we make a Ansatz for the 7 £ „ [which b y (|Cl | ) and 
( C2 ) is equivalent to the Bogoliubov Ansatz (|64[)]: 



lev = X B enu ~f nu (y = ±) 

n£{d,fc>0} 



(C7) 



It suffices for the Ansatz to be linear, since H x is 
quadratic in 7n„'s, and for it to be diagonal in the in- 
dex v, since both flC5| ) and ( p6[ ) are purely off- diagonal 
in Since the orthonormality conditions (|C4|) imply 



(" = ±), 



(C8) 



the -B Eru /s are orth ogon al matrices [with matrix indices 
(e,n)], so that Eq. ( |C7| ) can trivially be inverted: 



7w 



(C9) 



e>0 



We can deduce S Eq eve n wi thout having determined the 
B Env yet, by inserting (C9) into H x to transform ( |C5| ) 
into (C£): since both equations are off-diagonal in v, no 
diagonal terms levlev (= 1 ; i- e - n o constants) can arise, 
so that the constants in both equations must be equal; 
this yields Eq. for SE G . 



2. Determination of _B en „'s and e's 

To de term ine the coefficients B enl/ , we substitute the 
Ansatz (C7) it into the Heisenberg equation 



[%u,H x ] = vie {v = ±) 



(CIO) 



[which follows from (C4) and (C6) and is equivalent to 
(p5|)], evaluate the commutator using ( |C3| ) and (|C5j) , and 
equate the coefficients of j ni/ . This readily yields: 
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eB ek+ = kB £ - k _ + 2V- k B ed _ , 

sB e d+ = \Sd\Bed- > 

eB £d _ = \e d \B £d+ +J2w- k B e - k+ 

k>0 



(Clla) 
(Cllb) 
(Cllc) 
(Clld) 



We consider the e 7^ and e = solutions separately. 



a. Solutions 

For e / we write -B E£ 2- = q(e) £, where g(e) is 
a nor maliz ation factor to be determined below. Then 
Eqs. ( HI! ) yield Eqs. (|69|) afte r some simple algebra. 



Substituting these into (Clld) yields the eigenvalue 
equation 



S 1 {e)=e 2 /e 2 -l 



(C12) 



where the infinite sum Si(s) can be evaluated as follows: 

(C13) 



4F 2 



k 2 



k>0 

V - W //% (C14) 

(47rr/e)tan[7r(e/Ax,-P/2)]. (C15) 



Equating Eqs. ( |CT^ ) and ( |C15| ) gives the eigenvalue equa- 
tion (|6^) that determines the allowed e-values. 

The normalization factor g(e) can be det ermi ned by 
writing the e — e' 7^ 0, ^ = + version of Eq. (C8) in the 
form [using Eq. m£ 



(C16) 
(C17) 



e 2 {e) [e 2 d + e 2 S 2+ (e)] =1 
4V- 2 £ 2 



k>Q 



(k 2 -e 2 ) 2 



Noting from Eqs. flC13| ) and pL7| ) that 6*2+ (e) = 
ieSeiSi^), evaluating this by reexpressing the deriva- 
tive of — ( plq ) in terms of tanJJ. an d sim plifying using 
tan[] = (e^~P)/(7reT) [by (|5l5|)=(|C12|)1, one readily 
finds Eq. © for g(e). 

Note that the eigenvalue equation (^) is symmetri- 
cal under the transformation e — > —e an d the refore also 
has negative roots. However, from Eqs. ( |Cll[) the corre- 
sponding coefficients are given by B_ env = vB enu , thus 
the excitations corresponding to e and — e are not in- 
dependent, but related by <5_ e = a\ . This confirms 
that only non-negative eigenvalues need to be considered 
[as was intuitively obvious already when writing down 
Eq. (§)]. 



b. e = Solutions 

There are three situations in which a root e^p of the 
eigenvalue equation ( |67| ) can be equal to zero: 

(i) For arbitrary T, the lowest root in the P = sector, 
£0,0, approaches if and only if e d — ► 0. For e d 7^ 0, its 
asymptotic behavior close to or Al/2 is as follows: 



£0,0 ^ [£rf| 

A L ~ [A|+4^ 2 rA L ] 1 /2 

£M i _ 2TAl 
A L ^2 £ 2 -A|/4 



for 



£0,0 

Ar 







(C18a) 



for (C18b) 
Ai 2 



For e d = eo.o = 0, the coefficients B eo 0= o nv are simply 
the £d — > limits of the B Eo onv of Eqs. (|69|), 



<5 rad -(2Vyfc)^ 

[i + 4^ 2 r/A L ]i/2^ 



(C19) 



reflecting that 7^+ decouples from H x for e d = 0. 

(ii) The lowest root in the P = 1 sector, £0,1, identi- 
cally equals for all e d = and T. Solving the e = 0, 
r =/= versions of (Cll) directly (since (pa), derived for 



£ 7^ 0, can not be used here) gives 



B 



eo,i"+ 



[4y 2 + £ 2 ]V2 



(C20) 



reflecting that 7fc =0 _ decouples from 7?^. 

(iii) The second-lowest root in the P = 1 sector, £1,1, 
approaches zero if and only if T — > at £d = (to be 
precise, in this limit ei t i ~ 2v / TAi); this reflects the fact 
that H x has two zero modes for e d = T = 0, namely «o 
and a^. Note, though, that in this limit 7 El lV [found 
using ([39|)] and 7 eo do not reduce simply to 7^ and 
jk=o,u but to linear combinations of these (in Table §, 
this is indicated by braces). In the opposite limit of 
r/A L -> 00 at e d = 0, one has £14 ~ A L (i - ^fy). 



3. Consistency Checks 

Several consistency checks on the above solution are 
possible. Firstly, let us check Eq. (jCq): In the spe- 
cial case that e or e' is £o,i = , Eq . (C20) is easily 
checked to be consistent with Eq. (C8). If £ and e' 7^ 0, 
one finds by writing out Eq. (C8) for (e — e' ,v = — ) 
and (e 7^ e ', v — ±), respectively, that [analogously to 
T0),(C17)] the following relations must hold: 



1 



Q 2 (e)e 2 
r 2 



4V? fc 2 



(fc 2 -£ 2 ) 2 



fc>0 

4K 2 ££' 
k 



k>0 



E 

fc>0 



(fc 2 -£ 2 )(fc 2 - 

AV- 2 k 2 
(k 2 -e 2 ){k 2 - 



S a _(e)], (C21) 
^3+ (£,£')], 
S3- (£,£')]• 
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One can v er ify t hat indeed they do, by noting from 
Eqs. dgH), (|C17b that the sums defined above can be 
rewritten as 52- = S± + S2+ and 



S 3 ±(e,e') = 



1 



e5*i (e) - e'Si (e') eSi (e) + e'Si (e 7 ) 



and simplifying these using Eqs. flCl^ ) and jCl^ ), 

Secondly, one can verify explicitly that our transfor- 
mation does indeed diagonalize H x : insert the inverse 
Bogoliubov transformation (C9) and Eqs. ( |69| ) for the 
coefficients B env into the original form QC5| ) for H x and 
express the resulting y% sums in terms of Si, S2- and 
S3-: 

F a = M/2 + E fc / 2 + E^ 



fe>0 



fc>o 



E i %+ ;: fe'-Q{e)Q{e')e' 



e 2 Si(e) 



(C22) 



e,e'>0 

+ e 2 (s e£ ,S 2 _(e) + (1 - 6 BS ,)S 3 _ (e,e') 



(The terms with e or e' = can be checked to be zero.) 
Evaluating this using Eq. ( |C12 ) and the equations for 
Sg- and S3_ readily yields the sought-after diagonal form 
© for P K and confirms Eq. §§ for 5P G . 



4. Ground State Energy Shift 5E a 

We now show how to calculate the ground state en- 
ergy shift 5 E G ee 5E% + PSE^ of @, for T/A L > 1 
and both T ft = and T h /A L > 1, i.e. w e derive Eqs. @ 
and (|7§|). As explained in Section [VP , the coefficient of 



the P-dependent term, <5P<3, must be extracted with care 
to obtain the correct finit e-siz e spectrum. 

In the notation of Eq. (71a), Eq. (pq) becomes 



6E G 



E 



JV max 
j=0 



E- 5 



:,/| ■ --„./• ■ ■ -Al > . 'V./- (C23) 

where we introduced a "band cutoff" N max = D/Al to 
regularize the sum (with D ~ 1/a) . The task at hand is 
to perform the sum on j sufficiently carefully to extract 
its leading P-dependence. 

a. Zero magnetic field 

We first consider the case T / A L > 1 and T h — To 
isolate the P-dependence of 8j t p, we write Eq. ( 71b ) as 



Sj, P = 1/2 + g(j - P/2 + 6j, P ), 
A L (x-l/2) 



g(x) = arctan 

7T 



4tiT 



(C24) 
(C25) 



and solve (C24) for 5j t p by expanding its r.h.s. in the 
small parameter (<5j,p — P/2)Al/T <C 1, finding 



l/2 + g(j) P (A\ 



(C26) 



The first term is P-independent and gives the leading 
contribution to 6E G . The second term is O (Al/T), con- 
tains the full P-dependence of 6j t p to t his o rder and con- 
tributes to SEq. Inserting (C26) into ( C23| ) gives 



8E% = - 



A L N ™ 1/2 + g(j) 



(C27) 



= -r 
= -2r 



/ dy[n- 2 arctan y] + 0(A L /T) 
Jo 

ln(P/47rr) + I + O (Al/T, T/D) 



(C28) 



while the P-dependent part, SEq is equal to 

Nmax 

- £ o,p/2+(A l /4)£ 5 '(j) 

« -e ,p/2 + (A L /4) [.g(iV max ) - 5(1)] 
«-Ai[l/8 + 0(A£/r)] , 



(C29) 

(C30) 
(C31) 



where for (|C3l|) we used e ,p = for T ft = [by @]. 



6. Large magnetic field 

Next we consider the case T/Al 3> 1 and T^/A^ ^> 1 
(for arbitrary Th /T). This can be treated analog ously , 
except that now (C25) must be replaced by [from ([711])] 



, . 1 1 
g(x) = — arctan — 

7T 4-7T 



Ai(s-l/2) 



A L (a;-l/2) 

r 



[thus g'(x) is of order 0(A L /T, Al/Tu) for all a > 1 ], 
Since now e . P = A L (1 - P)/2 [by @ and dC18b|) ], 
Eq. ( |C3C| ) now yields A L [l/2 + 0(A L /r, A L /T ft )] , which 
is P- independent . This implies that SEq = to this 
order, in other words that for Th, T 3> A^ (and indepen- 
dent of the ratio Th/T) the ground state energy shift SEq 
is no longer P-dependent. 

The P-ind epen dent part of the shif t, 5E %, can be ob- 
tained from (1C271 ), plus the |e d |/2 of ( |C23|) : 
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5E% = M _ 



dy 



7r — 2 arctan y 



16n 2 Ty 



0(A L /T,A L /T h ) 



2T 



\n(D/\e d \) + l + 0(A L /T,A L /T h) T/D,T h /D,T/T h ) 



Although the integral can be evaluated for arbitrary values of its parameters using 



dy arctan(y — b/y) = yarctan(y — b/y) — — ln(y 4 + y 2 



we gave in (C32) only the limit of large magnetic fields, 
T h /T » 1. 



5. Fluctuations in M x 

The results for (N x ) and (A/ -2 ) discussed and inter- 
preted in Section VF are obtained as follows. We con- 
sider only e d = and the states |0)s oxt and o4 P |0)s ext , 
which represent, respectively, the physical ground states 
in sectors with excitation parity V ^ = and 1. We ex- 
press N x through aj. and /3g , using (^l|) and (|5^) : 

A = J2 ^tPk - Phk) + P0-/2 - ajao) . (C34) 

fc>0 

Now, in a P = 1 sector, we have [analogously to ( C43| )] 



1/2 - (ala ) 



1 
2 

0. 



i^eo.iO+^eo.iO- + E B e0+ B s0 ^ 

(C35) 



£>E0,1 



for all T/A L [using (|69|) and (|C20|)]. (Here and below 
the upper or lower signs in i (and = p) refer to P^, — 
(<5j p a eo p ) — or 1.) Since moreover /3fc|6) 5cxt = 0, we 

conclude from ( p34| ) and flC35| ) that (AQ =0 for both 
P = and 1 and all T/A L . 

The calculation of (A"" 2 ) is more involved: 



<A4 2 )+P/4 = 5>!a,> 



k>0 



2 E f 1 ^ B eo,pk+ B e ,pk- 



,k- 



k>0 



j>0 



2 Wfc* 



k>0 



e>0 



(e 2 - fc 2 ) 2 

2 



£>0 



fc>0 



(fc + e) 2 



fc>0 e>0 



For the first equation we u sed (C34), ( C35Q; for the sec- 
ond ©, @, @, dC||); for the thkd^), (|C19|) or 
( C20Q ; we tamed the "divergence" at e ~ fc using 

efc P 1 1 



(e 2 -*; 2 ) 2 (fc 2 -e 2 ) 2 2(fc 2 -e 2 ) 2(fc + e) 2 



2y 2 b- 



u2, 1 771 7 /l + 2y 2 -2fo\ 

b 1 ) V46 - 1 arctan = , 



V46~" 1 



(C32) 



(C33) 



and p erform i ng th e sum s over the first t wo te rms 
using (|C12|) , (|C2l| ) and ( |C13|) , thus obtaining ( |C3^ ), in 
which the last two (diverging) terms cancel exactly. 

The limit T/A L -» of (|C36|) yields (A/|) -» P /4 [fo r 
P = 1 the mode £1,1 — > 2vTA7 (compare Section |C 2 b| ) 
makes a non-zero contribution]. To obtain the leading 
behavior of ( p36 ) in the opposite limit T/Al 1, we 
evaluate the sums as integrals: 



de ■ 
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±£ 2 + 47r 2 r 2 



dk- 



2T 



(k + ef 



— 2 ln(TL) + 0[(A L /T) } . (C37) 



6. Perturbing around EK line by 5H' Z 

In this section we determine the scaling dimension 
jsx z of the operator 6H' Z of (|38|), which arises as soon 
as one leaves the EK line, i.e. when X z = 1 + 5\ z . To 
this end we perturbatively calculate, in a given subspace 
Syhys ( St , A~ c , A/V ) at = 0, the level shifts 5£{L) [see 
(P0|)1 induced by SH' Z . As first step, we express 5H' Z in 
terms of operators that diagonalize H'(SX Z = 0): 



d x ip s (x) = A L V T h i ( b is - b U) 

n q >0 

Sf s S z = S T S Z - 1/4 , 



(C38) 



(C39) 



S z 



B e , d _(4 + a e )(al,-a e ,) . (C40) 



[Eq. ( p38| ) follows from (|l4j| and (|2§ ), Eq. (|C39|) from 
0j, and Eq. (|C40| ) from (|f), |>8cf), and 
( C2\ ) .] These relations show that although the NFL fixed 
point spectrum is highly degenerate, there is no need for 
degenerate perturbation theory, because (E\8H' Z \E') = 
whenever \E) and \E') are degenerate but distinct eigen- 
states of H' [by inspection; compare Table |]. To first 
order in (5A 2 , the dimensionless energy shift 5£{L) of (|9C 
due to 5H' Z is thus simply given by 



S£(L)= (E\6H'\E) - (E m J5H'\E min ) /A 



6\ Z S T 



(E\S Z \E) - (E.\S z \E mh! 



(C41) 
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(Note t hat s ince d x (p s [by ( |C38| )] is linear in b qs and 
¥ gs , in ( |C4l| ) only in the second, "zero-mode term" of 
6H' Z contributes, which does not occur in the contin- 
uum limit considered in Ref. ^.) Now, i? m i n corre- 
sponds to the physical ground states in the two sectors 
S phys (S T = ±1/2, M c = 0,7V> = 0) [with £ NFL = in 
Table fif , namely 



|0)s cxt (l/2,0,0) , 5 i ,ol (3 )cS ox t(-l/2,0,0) 



(C42) 



From (|C40| ), one readily finds for these, respectively, 



S T (S Z 



£o,od+B £oo0 



j>0 



(C43) 



Evaluating this using ( |69a| ) and ( C19| ), the first term 
yields |[l + 47r 2 r/Ai] _ 1 1 for both cases, and the second 
term vanishes for Sd = 0. Thus the two- fold ground state 
degeneracy persists, in keeping with the fact that 5H' Z 
respects spin reversal symmetry. 

The first four excited states in Table |, with = g, all 
have St = 0, hence (StS z ) = 0. They are thus not 
shifted by 5H' Z themselves. Nevertheless, their relative 
shift w.r.t. the £nfl = physical ground states just dis- 
cussed is non -zero , since the ground states were shifted 
upwards: by (C41), it is 



S£(L) 



4vr 2 r/A L ]- 1/2 ~ L~ 1/2 . (C44) 



It follows [from (90)] that the sought-after scaling di- 
mension of 8H' Z is "fs\, = 1/2. Thus this perturbation is 
irrelevant, and the RG flow in the vicinity of the EK-line 
is always towards it. [It is easy to check that all 10 of 
the next-higher excited states in Table [i], with f nfl = \ ■, 
have S£(L) = -8\ z ±[l + 4 7 r 2 r/A L ]- 1 / 2 , which again is 
~ L -1 / 2 , as expected.] 

Let us now turn on a local magn etic field Sd = hi, in 
which case the second term of ( |C43| ) is non-zero and con- 
tributes to lifting the ground state degeneracy. In the 
continuum limit L — > oo (so that -C hi,T), it in fact 
ayes a much larger contribution than the first term of 



namely 8£(L) = ±^I(h), where h = \hi\/(2nT), 
and 1(h) is given, after the substitution x — e 2 /e|, by 
the following integral: 



1(h) = 



dx-- 



h 2 (x-l) 2 /2 + 2x 



(C45) 



h — >oo 



'-h\n(h/2) ■ 
2 

irh 



APPENDIX D: THE SINGLE-CHANNEL KONDO 
MODEL 

This Appendix deals with the anisotropic single- 
channcl Kondo (1CK) model, which is of interest not 
only as the most basic Kondo model, but also since it 
is equivalent to a dynamic two-state system coupled to 
an ohmic environment]!^! We shalLsaLye the 1CK model 
along the so-called Toulouse line,E3i£j the 1CK analog 
of the EK line, calculating the crossover of the finite- 
size spectrum from the free Fermi liquid fixed point to 
the strong-coupling Fermi liquid fpted point, well-known 
from Wilson's NRG calculations.E3 Since the 1CK cal- 
culation is a straightforward adaption of that developed 
above for the 2CK case, it will be presented in less de- 
tail than the latter, though technical differences will be 
pointed out. 



1. Conserved Quantum Numbers 

The 1CK model is defined by Eqs. (g) or @ for H , 
(||) for i?j n t and (j^) for Hh, the only difference being that 
the channel index only has the value j — 1 and hence can 
be dropped throughout. To exploit the fact that the to- 
tal charge is conserved, we transform from the a = (t, 4) 
basis to a y = (c, s) basis by writing: 



<Pc/ s (x) = (4>t(x) ± 4> l (x))/V2 . 
N c /s = (Ni±N l )/2, 



(Dla) 
(Dlb) 
(Die) 



Note that the normalizati on co nstants in Eqs. (Dla) and 
|Dlb| ) differ from that of (Die) [this contrasts with the 
2CK case, and affects many of the equations below]: the 
l/y/2 in the former ensures that the transformations for 
&| y and ifiy are unitary, so that these operators satisfy 
commutation relations analogous to those of b qa and 4> a 
[namely (U) and @]; the 1/2 in ([Did) ensures that N c 
and N s can be interpreted as half the total charge and 
the total electron spin, whose eigenvalues are are either 
both integers or both half-integers (whereas a \/\/2 in 
(Die) would have yielded irrational eigenvalues): 



M = (M C ,M S ) e (z + P/2) 2 



(D2) 



Here the parity index P equals or 1 if the to tal n umber 
of electrons is even or odd, respectively. Eq. (D2) is the 
free gluing condition for the 1CK model. 

Evidently, the total charge M c and the total spin, 

St = M s 



S z 



(D3) 



are conserved, where (D3) will be called the spin- 
conservation condition. Hence we can restrict our at- 
tention to the invariant subspace 



S phys (S T ,M c ) = {\M Cl S T -k'A)® WcSt + iA)} 



(D4) 
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The difference between Eqs. (D4) and ( p5| ) makes ex- 
plicit a major difference between the 1CK and 2CK mod- 
els: though for both the quantum number Af s fluctu- 
ates "mildly" between St T 1/2, the 1CK model lacks a 
"wildly" fluctuating quantum number such as Af x ; this is 
the "deep reason" why it also lacks NFL behavior. 

In the new charge-spin basis, the bosonizcd form of the 
1CK Hamiltonian takes the following simple form: 



H 


= A L [M c (l - P ) 


+ N 2 +J\ 










b ls b qs) , 




(D5) 












H z 




+ A L Af s 


S z , 


(D6) 


H x 


2a 


f^t + 


i.e. . 


(D7) 



S z [by (Dll)]. Thus it is trivial to find the spectrum, 
which turns out to coincide with the fixed-point spec- 
trum shown in the strong-coupling limit of Fig.^. Note, 
incidentally, that at the decoupling point the modei-can 
be mapped to a two-level system without dissipationBcil 
The other solvable point is the Toulouse point, with 



At s 2 - -Jl , 7 



* = V2-1 , 
^'W + h.c.) 



(D12) 
(D13) 



We henceforth focus on this point, which is the analog of 
the EK line in the 2CK context, since the factors e %Vs can 
be treated by refermionization, as shown below. Note, 
though, that the spin —j*S z /V2 from the conduction 
band that is tied to the impurity does not fully screen 
the latter. 



2. EK transformation 

To simplify H z , we use the same Emery-Kivelson trans- 
formation U(\) = e* s *V'W of (j3TJ) as for the 2CK model, 
but now with (p s (x) given by Eq. ( jDlb| ) instead of @. 
The impurity spin operators S+ and the spin field ip s 
transform according to Eqs. (|32|) and (|34|), just as in the 
2CK case, but in contrast to the latter, 



U(H + H Z )U- L =H a + (A 2 /V2 - 7)^^.(0)5, 
+\ z A L Af s S z + const , (D8) 
ip a (x) -» ip a (x)e iV2a ^ s ' "«*»('»/«) , (\ x \ < L) . (D9) 

Moreover, since in the 1CK case the spin density is 
d x ip s (x)/(2Tr^/2), a spin —~/S z /V2 from the conduction 
band is tied to the impurity. To eliminate the S z d x tp s 
term, we choose 7 = \ z j\f2 (in contrast to 7 = X z for 
the 2CK case). Then 

H'(X ± = 0) = A L [ti c {N c + 1 - Pa) +M 2 S + \ Z N S S Z 
+ « ( h \c V + + H h + const , (D10) 

q>0 

and H' ± contains the factors e ±*(V2-A*A/2)<^(0) _ 

These factors simplify for two special values of X z . The 
first case, A z = 2 (i.e. 7 = V2), is called the decoupling 
point, since the (^-dependence drops out completely: 



H' x = ^(S + FlF t +h.c.) 



2 a 



(Dll) 



In this case, the spin —^S z /V2 from the conduction band 
that is tied to the impurity is precisely —S z , thus we have 
perfect screening. Indeed, by ( p9| ) the phase shift S in 
V-'qCO - ) = e l2S ip a (0 + ) is \S\ = tt/2, corresponding to the 
unitarity limit. The dynamics of the electron- hole excita- 
tions described by the c/? SjC fields evidently decouples from 



3. Refermionization 

To ensure proper anticommutation relations for the 
pseudofermions to be defined below, it is convenient to 
make one more unitary transformation with the opera- 
tor U2 = e l7rA ^ s * , which changes the phases of the Klein 
factors and the spin operators (and of the basis states in 
Fock space): 



u 2 fIf^ 



-" s -fIf ] , 



s± 



(D14) 
(D15) 



Then H' x of Eq. (|>T|) takes the very simple form 

U 2 H' A U^ = ±= (4V,(0) + iPt(0)c d ) , (D16) 



where we introduced the following pseudofermions, 

4 = S + e^°~ s ^ , 4c - S z + 1/2 , 



ip s (x) 



I T ~i(Ar s -s f rn(S T )\l+P]/2)'2nx/L-iv s (x) 



a 1 / 2 



(D17) 

(D18) 
(D19) 



with sgn(Sr = 0) = 1 . By i nclu ding the factor 



in the definition (D18) of ip Sl we pur 



s i[l+sga(S T )P]ltx/L 

posefully ensured that ip s has the same boundary condi- 
tions (namely periodic) for both P = and 1, in order 
not to have to distinguish between these two sectors (the 
reason for the sgn(Sr) facto r is e xplained below). As a 
consequence, the k's in Eq. ( D19 ) must be of the form 



k = A L [n- k -sgn(S T )l/2], n- k e Z. 



(D20) 



[i.e. the periodicity parameter P of (Q) here equals 1]. 
By Eqs. (p9h, these pseudofermions have the properties 



2G 



{ c fcs' C \' S } — ^kk' > 

{cd,4) = 1 > 

[c d ,A/" s ] = c d ■ 



(D21) 
(D22) 
(D23) 
(D24) 



Note in particular that the anticommutation of a and cj, s 
is ensured by the factor e~ %7!l ^ s in the definition ( |D 1 7| ) of 
Cd- Note further that the individ ual action of both Cd 
and ip 8 violates the conservation ( |D3| ) of the total spin 
St = M s + S z . When diagonalizing H' , we shall therefore 
work not in the physical subspace 5 p h ys (5'T, M c ) of (D4), 
but in a correspondingly extended subspace S cx t(M c ), in 
which M s is unrestricted and not linked to S z . At the 
end of the calculation we shall retain only the physical 
states in S ex t (yV c ) , which we identify using a generalized 
spin-conservation condition to be derived from ( |D3| ). 

Let |0),s be a free reference ground state in S ext de- 
fined as in (pQ), and let : : denote normal ordering w.r.t. 

it. Then : c\cd ■— 4 c d — rr d °'. Moreover N s = J2k '■ 
which counts the number of s-pseudofermions, is 



c\ c 

ks 



ks ' 



related to M s by [compare (51)] 



N s =M S - sgn(S T )P/2 



(D25) 



The sgn(Sr) factors above are needed because we pur- 
posefully included one in the refermionization relation 
(|D18j) ; we did this to ensure that the spin reversal trans- 
formation (S Z ,M S , ifs) — > (—S z ,—M s ,—<p s ) can also be 
simply implemented in terms of the new pseudo-fermions, 
for which it implies 

(S T ,c d ^s,Ns) -> (S T ,c d ,iPi-N s ). (D26) 
The pseudofermions' kinetic energy is [compare (|5S 



y 

k 



/ , k :C ks C ks : 



&LNJ2 



E 

Q 



q b\ s b qs . (D27) 



This result can be used to rewrite Hq in (D5) in term s 
of the new pseudofermions Cj, s and Cd- Though (D27) 
differs from ( |D5| ) by terms in both A/ 2 and M s , the dif- 
ference can be expressed in terms of c\cd using the spin- 



c\c d 



conservation condition (D3), namely M s = St + ^ 
(those states in 6> cxt for which (|D^) does not hold will 
be discarded at the end anyway). In this way, the EK- 
transformed Hamiltonian of Eqs. dD5| ), (D8) and (D16) 
can be brought into the following refcrmionized form: 



U 2 H'U^ 1 = H c + H s + E G + const. 



H c = J2q b\ c b q 
<? 

H s = ^ k : c\ s c- ks : + e d : c\c d : 

k 

+ v/A^f^(ci s c d + 4cfe s ) 



(D28) 
(D29) 



(D30) 



Er 



M c (M c + l-P ) + ±S T [S T + sgn(S T )P] 



i-p 

8 



4 



Sxh e 

fir , 



Ed 



n (0) - i 

n d 2 



£d — ^-L^dfi + hi 
£ d ,o = (K - 1)St - sgn(S T )P/2 . 



(D31) 

(D32) 
(D33) 



Evidently, the charge sector decouples completely. In 
the spin sector, H s corresponds to a quadratic resonant 
level model, whose "resonant level" has energy Sd and 
width T = A^/4a, and Eg is the "free ground state en- 
ergy" of the spin sector in the extended subspace <S ext , 
in the presence of magnetic fields. Note that in the ab- 
sence of magnetic fields, i.e. hi = h e — 0, we also have 
sgn(e ( ;) = — sgn(SV); thus it follows by inspection that in 
this case H s and E q are invariant under the spin rever- 
sal transformation ( D26| ), as they should be. As for the 
2CK case, we henceforth set h e — 0, the generalization 
to h e ^ being straightforward. 



4. Diagonalization of H s 



We wish to bring H s into the diagonal form 



if, 



E 



e :clc F 



5Er 



(D34) 



Here : : denotes normal ordering of the c £ operators with 
respect to an exact reference ground state |0)s cxt of the 
subspace <5> GX t, defined by the conditions c e |0),s cxt = for 
e > and c^|0)5 oxt = for e < 0. This may be accom- 
plished by a unitary transformation of the form 



= ]T B En cl. 



(D35) 



n—{k,d} 



Analogously to Appendix ^| the -B en 's and e's can be 
determined starting from the relations 



[H s ,cl] = £ 4, {C e ,4'} = <w 
One readily finds the following results: 



£-£d 
B 



ed = 



COt(7T£/A L ) , 



B 
SE 



r 2 7r 2 +rA L + (e-e d ) 2 
f 



1/2 



B 



ed 



E £_ Y. k ~ £dn d 



(0) 



(D36) 

(D37) 

(D38) 

(D39) 
(D40) 



£<0 



fe<0 



The eigenenergies e are the roots of Eq. ( D37 ). Their 
general behavior as functions of T and Ed can be deter- 
mined graphically, similarly to Fig. |^. To identify the 
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crossover scales of the problem, we write a general solu- 
tion, in analogy to (\n\), as 



sj = A L (j- 1/2 + Sj) 



5j = — arctan 

7T 



Trr 



(D41a) 
(D41b) 



Evidently, all solutions with \ej — ed\ <^.Y are strongly 
perturbed, with upward or downward shifts Sj ~ ±1/2 
for £j > or < e^, whereas those with \sj ~ Ed\ ^S> Y are 
only weakly perturbed, with Sj ~ 0. (The solution Sj 
closest to Ed can be associated with the d level, which, 
measured in units of A^, is pushed to the integer closest 
to it as F/Ai — > oo.) 

This implies the following crossover scales: 

(i) Without magnetic fields (hi = h e = 0, so that Ed — 
A-L^dfi) and in the limit Y ^> Al (i.e. also Y 3> |£<i|), 
the crossover scale separating the strongly and weakly 
perturbed spectral regimes is T, which can thus again be 
associated with the Kondo temperature, i.e. Tk — Y. 

(ii) For a large local magnetic field hi S> A^ (so that 
£d — hi), the /li-induccd shifts in the lowest-lying levels 
with \sj\ <C r become large (~ 1/2) roughly when hi 
reaches the crossover field h c ~ Y. In other words, for 
hi ^> h c , the local magnetic field is strong enough to ef- 
fectively erase the effects of spin-flip scattering from the 
lowest-lying part of spectrum. Note the contrast to the 
2CK case, where the crossover field is smaller, namely 
h c ~ VTAT. 

For r 3> ed,A.L, the ground st ate energy shift SEq, 
calculated similarly to Appendix C 4 , turns out to be 
SE G = 5E G + SE G , where 8E% w -Yln(D/Y) is the 
"binding energy" , and 



SE G 



£d[l/2- 



,(o) 



(D42) 



where the dots represent terms that are either indepen- 
dent of P and £ d , or of order 0(e 2 JY, A 2 L /Y, A L e d /Y). 
Note that SE G precisely cancels the Ed term in Eq of 
( |D3l| ), thus the fixed point spectrum at Y/Al = oo sat- 
isfyingly does not depend on the parameter A* = 2 — \/2 
occurring in Ed- 



5. Generalized Spin-Conservation Conditions 

To identify and discard all states in 5 ex t that violate 
the total-spin-conservation condition (D3), we now derive 
a generalized spin- conservation condition. The argument 
is analogous to that for the 2CK generalized gluing condi- 
tion in Section V C , but more straightforward, since H s 
conserves the number of c n excitations (not only their 
number parity). 

The number of excitations of a general eigenstate \E) 
of H s relative to 5 ex t is 



WE 



ClCr 



\E). 



(D43) 



When r is turned off adiabatically and \E) reduces to 
\E) = limr— >o its excitation number Af§ remains 
fixed. It hence equals ^^(r — > 0), which can be writ- 
ten as 

(E\Y^:c\c- ks : + :c\c d : \E) 

k 

- (E\ [N s + S z + 1/2 - t4 0) ] \E) . 

Using ( |D25| ) for N s and imposing the condition that any 
physica l \E ) must satisfy the total-spin-conservation con- 
dition (pa), we obtain 



E 



S T ~ [sgn(S T )P - l]/2 
St- [sgn(S T )P + l]/2 



This generalized spin-conservation condition specifies 
which of all the possible states in S P y .t a re physical; it 



supplements the free gluing condition (D2), which stipu- 
lates that St T 1/2 must be integer (half- integer) if JV C 
is integer (half-integer). 



6. Finite-Size Spectrum 

We consider here only the case Pq = 1 of anti-periodic 
boundary conditions (Pq = is analogous), and zero 
magnetic fields, hi = h e = 0. The construction of the 
finite-size spec trum is entirely analogous to the 2CK case 
of Section VE, but a little more cumbersome, since £d,o 
of (D33) and hence also Ed is not equal to zero; instead it 
depends on St, i-e. changes from one sector S ex t to the 
next. The results are summarized in Fig. [?] and Table III. 
The latter's caption also summarizes the technical details 
of the construction. 

(i) Phase-Shifted Spectrum: — The evolution of the 
phase-shifted spectrum £ p hasc for A z € [0, A*] at Aj_ = 
is given by H'(X± = 0) of ( p!0| ); it evolves linearly 
with increasing A z , from £f ree at X z = to £ p hase at 
A 2 = A* = 2 — y/2, as shown in Fig. @(a). 

(ii) Crossover Spectrum: — The crossover spectrum as 
function of Y/Al & [0, oo] at the Toulouse point A z = A* 
is shown in Fig. ^(b). The spectrum evolves continu- 
ously from the phase-shifted values £ p hase at T = to 
a Fermi liquid fixed-point spectrum £fl at Y/Al = oo, 
which is constructed analytically in Table [II. The fixed- 
point spectrum corresponds pr ecis ely to the Fermi-liquid 
spectrum of free fermions [of ( Pa )] obeying the periodic 
boundary condition Pq = 0. This agrees with the stan- 
dard results p£ Wilson's numerical renormalization group 
calculations ,Ej and is expected, because at Y — oo one 
electron is bound so tightly to the impurity that the to- 
tal number of free electrons effectively changes tw one, 
and hence the chemical potential shifts by Ai/2£3 

The fact that all £fl's are integers or half-integers is a 
very direct sign of Fermi liquid physics, since it implies 
the absence of non-fermionic operators. 
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Finally, it is instructive to deduce another well-known 
fact, namely that a local magnetic field hi is a marginal 
perturbation, from the deviations from £fl which it pro- 
du ces: F or the lowest-lying levels, the value of the shift Sj 
of (D41) at the Fermi liquid fixed point (hi = 0, T > A L ) 
is |(<5j)fl| = 1/2. For a small local field /ii < Ai, an ex- 
pansion in powers of hi/ Al and Al /T yields 



(Sj)fl 



hi_ 
7r 2 r 



(j - 1/2 + - £ , d ) 2 Al 



7r 2 r 2 



Since the L-dependence of the leading term is ~ hiL°, 
the local magnetic field has scaling dimension 7^ = [cf. 
(|90|)] and hence is a marginal perturbation (in contrast 
to the 2 CK ca se, where it is relevant with jf li = —1/2, cf. 
Section VI C ) . A marginal perturbation always implies 
the existence of a line of fixed points, parameterized by a 
non-universal quantity. Indeed, for hi non-zero, the fixed 
point spectrum obtained for L — *■ 00 is non-universal, 
since in this limit ( D41b| ) yields non-universal shifts for 



the lowest- lying levels, namely Sj — > — — arctan(-7rr//i 
In contrast, for the 2CK case with T, hi 7^ 0, the limi t 
L — > 00 necessarily implies T,T/j ^> Al; hence ( |71b| ) 
yields universal shifts Sj^p — > 1 [see Fig. Q], which is why 
along the EK line the phase-shifted fixed point spectrum 
£ p h of Fig. ||(c) is independent of T and hi too (though 
away from the E K line it does acquire a a dependence on 
hi/T, see Section VIC). 
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(N c .N s ,N f ,N x ) = (0,0,0,0) 



(0,1,0,1) 



(0,0,0,2) 



FIG. 1. Under a succession of spin flips, A4 fluctuates 
mildly between St T 1/2 (here St = 1/2); in contrast, M x 
fluctuates wildly, since it can acquire any value consistent 
with the gluing conditions (^). The dotted line represents 
the reference energy up to which the free Fermi sea is filled 
for Po = 1, the filled and empty circles represent filled and 
empty single-particle states with energy k. 
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(a) P = 0, £ d =0 
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FIG. 2. The left and middle columns show, for var- 
ious combinations of (AT X ,S Z ), some representative free 
many-body states [eigenstates of H'(X± = 0) of (^)] in a 
general physical subspace S p h yB (ST, N C1 A//), constructed in 
terms of both c n 's and a n 's, thus illustrating the transfor- 
mation (^) between these operators. The braces between 
these two columns indicate that the states to their left and 
right are not in one-to-one correspondence, but linear combi- 
nations of each other. The right column shows some of the 
exact physical many-body eigenstates \E) of the full H' of 
(^) at T/Al >• 1. Each \E) is labeled by the excitation en- 
ergies Ej t p and kp of its occupied single-particle states. When 
r is turned off to 0, each state in the right column reduces to 
the free state in the same row in the middle column, unless 
they are separated by braces, in which case it reduces to a 
linear combination of the two degenerate free states grouped 
within the braces to its left. The first four rows in (a) and 
the two rows in (b) correspond, in that order, to the first four 
St = 1/2 states and the two St = states listed in Table |; 
the fifth row in (a) is the spin-reversed partner of the first row 
in (a), illustrating how the two- fold degeneracy guaranteed by 
spin-reversal symmetry comes about due to the presence or 
absence of a eo,p = excitation; (c) and (d) illustrate the 
case Ed 7^ relevant for non-zero magnetic fields. 
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FIG. 3. Graphical solution of the eigenvalue equation ( p7| ) 
for (a) P = and (b) P = 1. The vertical solid line marks the 
position of I Ed I . Dashed and solid lines represent the left- and 
right-hand sides of (|6t|) ; their intersections give the allowed 
eigenvalues e. The "amplitude" of the dashed lines is pro- 
portional to T/Al; if this increases from 0, the e's thus shift 
away from their free values \Ed\ or k = Ai,[nj. — (1 — P)/2]. 
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FIG. 4. (a) Evolution of the excitation energies Sj,p, found 
by numerically solving the eigenvalue equation ( |fj7| ) (or by a 
graphical analysis as in Fig. |3|). On the left the evolution is 
shown as function of F/Al S [0, oo) at Tj, = 0, and on the 
right as function of T^/Ax G [0,oo) at fixed r/Ax ^> 1, for 
(a) P = and (b) P = 1. These excitation energies are com- 
bined in Table | with excitations in the charge, spin and flavor 
sectors to obtain the evolution of the full finite-size spectrum 
shown in Fig. |(| 
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FIG. 6. Evolution of the many-body finite-size spec- 
trum of the 2CK model, for antiperiodic boundary condi- 
tions (Po = 1), from the free Fermi liquid fixed point to 
the NFL fixed point, and the additional crossover induced 
by a local magnetic field to a phase-shifted Fermi liquid fixed 
point. All eigenstates of H' of Eq. (JH^) are shown for which 
£nfl < 1 , as well as some higher-Jjjing states, with de- 
generacies given in brackets (in Ref.E3, the degenaries for 
£nfl = 1 were incorrect), (a) When A z is tuned from to 
its Emery-Kivelson value X z — 1, with Ax = £d = 0, the free 
Fermi-liquid spectrum £{ TCC at \ z — evolves smoothly into 
a simple phase-shifted spectrum £ phase at A z = 1. (b) When 
r/A_L = Ax/ (4a Ax) is tuned from to oo along the EK line, 
i.e. with X z = 1 and Ed = 0, the spectrum crosses over from 
£ P hasG to the non- Fermi liquid spectrum £nfl at F/Ax = oo, 
which agrees with NRG and CFT results, (c) Turning on a 
local magnetic field Ed = hi (with h e — 0) by tuning |ed|/r 
from to oo with A, = 1, T > Ax fixed, then induces a 
further crossover from £nfl to £ p h- For the lowest levels this 
crossover occurs when led | /T > 1, since then the crossover 
parameter used in Fig. |4j, namely 7), /Ax = (£d/F) 2 (T / Al), 
is ^> 1. The £ p h spectrum is identical to the phase-shifted 
spectrum £ p h aS e of A z = 1 and Ax = Ed — 0, apart from 
a degeneracy factor of two due to the lack of spin reversal 
symmetry. 
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FL 



FIG. 5. Sketch of the different Fermi liquid and non-Fermi 
liquid regions for a finite magnetic field on the EK line. 
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FIG. 7. Evolution of the many-body finite-size spec- 
trum of the 1CK model, for anti-periodic boundary condi- 
tions (Po = 1), from the free Fermi liquid fixed point to the 
strong-coup ling Fermi liquid fixed point. All eigenstates of 
H' of Eq. (D28) are shown for which £fl < 1, as well as 
some higher-lying states, with degeneracies given in brack- 
ets, (a) When A z is tuned from to its Toulouse-point value 
A* = 2 — v2j with Ax = £<i = 0, the free Fermi-liquid 
spectrum £{ Iee at A 2 = evolves smoothly into a sim- 
ple phase-shifted spectrum £ phase at A z = A*, (b) When 
r/A_L = Ax/(4oA_l) is tuned from to oo at the Toulouse 
point, i.e. with A z = A* and Ed = 0, the spectrum crosses over 
from £phase to the strong-coupling Fermi liquid spectrum £fl 
at T/Al = oo. The latter is identical to the free Fermi-liquid 
spectrum (X z = Aj_ = £d = 0) for periodic boundary condi- 
tions (Po = 0), in agreement with Wilson's NRG results. 
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TABLE I. Construction of the 2CK model's finite-size spectrum for Po = 1, corresponding to Fig. |(| The table shows all 
states that have excitation parity = [see ([75])] and a NFL fixed-point energy £nfl that is < 1, as well as some higher-lying 
states. (The states with P^ = 1 double the degeneracies of those with P^ — listed here, as explained below.) All energies are 
given in units of Al, e.g. £g = Eg/Al, with degeneracies in brackets. States in the same sector S p hy S (ST,Afc,Aff) are grouped 
together between a pair of horizontal lines and have the same Eg and SEq. (i) The construction of the phase-shifted spectrum 
for Ax = and Ed — is shown to the left of the brace column: in each sector, we list the lowest-lying free eigenstates 
of H'(\± = 0), some of which are illustrated in the left and middle columns of Fig. [| Each such state is labeled by the 
further quantum numbers (Af 3 ,Af x , S z ), satisfies the free gluing conditions (^), and has energy £f ICC or £ p hasc for \ z = or 1, 
respectively, as shown in Fig. |^(a). (N y = 0' here denotes a particle-hole excitation with M y = 0, {b\ y bi y ) = 1 and dimensionless 
energy q y — 1.) (ii) The construction of the crossover spectrum for Ai / at A z = 1 and Ed = is shown to the right of 
the brace column: in each sector, we list the lowest-lying physical eigenstates \E) of the full H'(X± 7^ 0), some of which are 
illustrated in the right column of Fig. ^| Each such j_E) is characterized by the excitation energies E ex = £j,p. kp or q y of those 
excitations ct\. p , (3^ or b\ y which it contains relative to the reference state |0},s ext in S esc t(ST,J^'c,J^'f) [see (|63J)], and satisfies 
the generalized gluing condition (f|). As r/A L increases from to 00, the excitation energies evolve from £ ex (Q) — > £ ex {oo) 
[as can be read off from Fig. Hfa)]; correspondingly, the energy of each eigenstate |-E) evolves from £ phase = £g + £ex(0) to 
£nfl = Eg £ex(oo) ± S£ G (the sum goes over all excitations listed), as shown in Fig. [](b). When T is turned off, each \E) 
reduces to the free state on its left in the same row, unless they are separated by braces, in which case \E) reduces to a linea r 
combination of the two degenerate free states grouped within the braces to its left (as illustrated in Fig. Hand Section C2b). 
By spin reversal symmetry, for each state shown here (all have = 0) there exists a degenerate partner with = 1, 
obtained by setting (S T ,JVf) — * {-St, —Mj) and either, for A± = 0, (N S ,N X ,S Z ) — > (—N 3 ,-N x ,—S z ), or, for X± / 0, by 
adding (or removing) a £0 p — excitation if it was absent (or present), as illustrated by the first and fifth rows of Fig. ^|(a). 
The degeneracies in Fig. |j(a) and (b), and summarized in Table |n| are thus twice those listed here, (iii) A crossover induced 
by a local magnetic field £d = hi (with h e — 0) occurs as Th/Al = Ej/(FAl) increases from to 00 (at fixed T/Al 3> 1), as 
shown in Fig. ^(c). It results from the further evolution of the listed excitations' energies £ ex with Th/Al [as can be read off 
from Fig. ^(b)]. The fixed-point spectrum at Th/Al = 00 is given by £ p h = £g + Y1 £ex(F/ 'Al = 00, Th/Al = 00) (note from 
(Q) that 5££ = for Th/Al > 1). 
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TABLE II. Summary of the finite-size spectrum of Fig. ^| 
for the 2CK model, at the four points A z = Ax = £d — 

(ffree); A* = 1, Ax = Ed = (fphase); A z = 1, T/A L = OO, 

e d = (£ N fl); and A z = 1, T/A L = oo, T h /A L = oo (£ ph ). 
We list all energies £ < 1 (in units of Az,) and give their total 
degeneracies in brackets. 
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TABLE III. Construction of the 1CK model's finite-size spectrum for hi = h e = and Po = 1, corresponding to Fig. [?]. 
The table shows all states with St > that have a strong-coupling Fermi liquid fixed-point energy £fl that is < 1, as well 
as some higher-lying states. (The states with St < double the degeneracies of those with St > listed here, as explained 
below.) All energies are given in units of Az,, e.g. £q = Eg/Al, with degeneracies in brackets. States in the same sector 
Sphys{ST, A/c) are grouped together between a pair of horizontal lines (and have the same Sd, Afg,, £ G and S£ G ). (i) The 
construction of the phase-shifted spectrum for Ax = is shown to the left of the braces: in each sector, we list the lowest-l ying 
free eigenstates of H'(X± — 0); each is labeled by the further quantum numbers (Af B , S z ), satisfies the free gluing condition (D2) 
and spin-conservation condition (D3), and has energy £f Tee or £ phase for A z = or A z = A* = 2 — v2, respectively. (JV V = 0' and 
the braces here have the same meaning as in Table |.) (ii) The construction of the crossover spectrum for Ax 7^ at A z = A z 
is shown to the right of the braces: in each sector, we list the lowest-lying physical eigenstates \E) of the full H'(X± 7^ 0) 
in that sector; each such \E) is characterized by the excitation energies £ ex = e, |e| or q y of the excitations c^ >0 , 5 £ <o or 
b\ y which it co ntains relative to the reference state |6)s cxt in S ex t(ST, M c ), and satisfies the generalized spin-conservation 
condition (D44). (For e < 0, [s| denotes the hole excitation c E |0)s oxt .) As T/Al increases from to 00, the excitation energies 
evolve from £ ex (0) — £ ex (oo); correspondingly, the energy of each eigenstate \E) evolves from £ p hasc = £g + ~}2£ex(0) to 
£fl = £g + ^2,£ex{po) ± S£q — 1 ~ 4 A (the latter constant corresponds to subtracting £ m in, the sum goes over all excitations 
listed), as shown in Fig. [?]. By spin reversal symmetry, each St > state shown here has a degenerate partner with total 
spin —St, obtained, for T = 0, by setting (Af s ,S z ) — (— Af s , —S z ), or, for T / 0, by setting (ed,J^^) — (— £d, — A/g) and 
interchanging particle- and hole-excitations, ej — e [cf. (D26)]. For all St 7^ levels, the degeneracies in Fig. m are thus twice 
those listed here. 
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